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EXACT  PERFORMANCE  OF  GENERAL  SECOND-ORDER 
PROCESSORS  FOR  GAUSSIAN  INPUTS 

INTRODUCTION 

The  performance  of  weighted  energy  detectors  and  correlators  for 
processing  deterministic  and/or  random  signals  in  the  presence  of 
nonstationary  noise  is  a  topic  of  frequent  interest.  Most  often,  a  second- 
moment  approach  is  adopted,  whereby  the  means  and  variances  of  the  decision 
variable  under  the  various  hypotheses  are  evaluated  and  employed  in  a  central 
limit  assumption  to  get  approximate  false  alarm  and/or  detection  probabilities. 
This  approach  is  suspect  for  small  false  alarm  probabilities  or  for  cases 
where  the  decision  variable  is  not  the  sum  of  a  large  number  of  independent 
random  variables  all  of  comparable  variance. 

A  recent  technical  report  [1]  has  presented  an  accurate  and  efficient 
method  for  evaluating  cumulative  and  exceedance  distribution  functions 
directly  from  characteristic  functions.  This  approach  is  very  fruitful  for 
determining  the  performance  of  general  time-varying  second-order  processors 
with  nonstationary  nonzero  mean  Gaussian  inputs,  since  the  characteristic 
function  of  the  decision  variable  can  be  evaluated  in  closed  form  in  these 
cases. 

We  will  consider  three  classes  of  processors  and  derive  the 
characteristic  functions  for  all  three  decision  variables  in  closed  form.  The 
first  two  classes  are  special  cases  of  the  third,  but  are  of  interest  in  their 
own  right,  since  they  include  and  immediately  reduce  to  many  practical 
processors  in  current  use.  Also  there  is  no  need  to  solve  for  the  eigenvalues 
and  eigenvectors  of  a  general  symmetric  matrix  that  is  encountered  in  the 
third  more-general  class  of  processors.  Rather,  the  characteristic  functions 
are  given  directly  in  terms  of  specified  processor  weights  and  input 
statistics. 

There  has  been  considerable  effort  on  this  problem  in  the  past;  for 
example,  see  [2,3]  and  the  references  listed  therein.  Most  of  the  lengthy 
analytical  derivations  and  results  have  been  aimed  at  getting  workable 
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expressions  for  the  probability  density  function  and/or  cumulative 
distribution  function.  Here,  when  we  consider  our  three  classes  of 
processors,  we  encounter  characteristic  functions  which  are  more  general  than 
that  given  in  the  recent  paper  for  a  filtered  analog  processor  [3,  eq.  5]; 
thus  specialization  of  our  results  will  yield  those  of  [3]  and  the  references 
listed  therein.  The  technique  employed  here  to  proceed  directly  to  the 
cumulative  and  exceedance  distribution  functions  is  a  numerical  one,  as  given 
in  [1],  and  does  not  require  any  series  expansions  or  analytical  manipulations 
at  all.  The  asymptotic  behaviors  of  the  cumulative  and  exceedance 
distribution  functions  on  both  tails  are  easily  observed  and  will  be  found  to 
corroborate  the  comment  made  in  [3,  p.  673]  that  these  tails  are  generally 
exponential  rather  than  Gaussian;  however,  there  can  be  a  considerable 
transition  region. 

The  programs  listed  in  the  appendices  require  the  user  merely  to  input 
his  processor  weights,  signal  constants,  and  noise  statistical  parameters  in  a 
series  of  data  statements  at  the  top  of  the  program,  and  to  select  values  for 

L,  limit  on  integral  of  characteristic  function, 

a,  sampling  increment  on  characteristic  function, 

b,  additive  constant,  to  guarantee  a  positive  random  variable,  and 

Mf,  size  of  FFT  and  storage  required. 

Selection  of  L  and  a  is  largely  a  matter  of  trial  and  error  and  is  amply 
documented  in  the  examples  in  [1]. 
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A  PARTICULAR  SECOND-ORDER  PROCESSOR 


Before  we  embark  on  the  analysis  of  the  particular  second-order  processor 
of  interest  in  this  section,  we  solve  the  following  simpler  statistical 
problem.  Let  s  and  t  be  real  jointly-Gaussian  random  variables  with  means 
ms,  m^.,  standard  deviations  as,  and  correlation  coefficient  p; 
thus  s  and  t  are  statistically  dependent.  Consider  the  random  variable 

x  =  as2  +  bt2  +  cst  +  ds  +  et,  (1) 

where  weightings  a,  b,  c,  d,  e  are  arbitrary  real  constants.  The 
characteristic  function  of  random  variable  x  is  defined  by* 


f  (?)  =  exp(ifx)  =  exp(if (as2  +  bt2  +  cst  +  ds  +  et))  = 


=  J|du  dv  exp(iif(au2  +  bv2  +  cuv  +  du  +  ev))  Pst(u,v),  (2) 

where  the  joint  probability  density  function  of  s  and  t  is 


Pst(u,v)  =  ^2ir  VtlfT7)  exp 


r  'u-msV  &  /v-mt\2  _  2pM\  /  v-mt 


“1 


2( 1-p) 


(3) 


Substitution  of  (3)  in  (2)  and  use  of  the  double  integral 


lldx  dy  exp[-ax2  -  sy2  +  2Yxy  +  ux  +  vy]  =  exp 


2  2 

Bu  +  av  +  2yu v 


4(aB-Y  ) 


for  ar  >0,  Br  >  0,  ar8r  >  Yr  , 


(4) 


*  Integrals  without  limits  are  over  (-«>,  +■<*>) . 
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(where  sub  r  denotes  the  real  part  of  complex  constants  a,  6,  y,  u,  v) 
yields,  after  an  extensive  amount  of  manipulations,  the  characteristic 
function  of  random  variable  x  as  the  compact  closed  form  expression 


f 


x 


exp 


if 


Nq  -  iTN1  -  ?2N2 
1  -  i?D1  -  $2D2 


(5) 


The  required  real  constants  in  (5)  are  given  directly  in  terms  of  the 
processor  weights  and  statistical  parameters  as 


D1  =  2(a a2  +  ba2  +  Cp0$0t)  , 

D2  =  (4ab  -  c2)(l  -  p2)a2a2  , 

2  2 

Nq  =  am$  +  bmt  +  cm$mt  +  dms  +  emt  , 

Nj  =  (4ab  -  c2)(-|  m2a2  +  j  m2a2  -  pm^a^)  + 

+  (2ae  -  cd)as(mtas  -  pmsat)  + 

+  (2bd  -  ce)at(msat  -  pmtas)  - 


,1  ,2  2  .  1  2  2  ,  .  , 

~  "2  d  °s  2"  6  CTt  dePasat) 


N2  =  -  (ae^  +  bd2  -  cde) (l-p2)a2a2 


For  later  reference,  the  mean  and  variance  of  x  follow  from  (5),  upon 
expansion  of  JJn  fxCf)  in  a  power  series  in  T,  as 


(6) 


ux  ~  No  +  2  Dr 


2 

1 


2N0Di 


(7) 
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(When  =  0  in  (6),  it  can  be  shown  that  Dg  <  0;  thus  characteristic 
function  (5)  never  possesses  any  singularities  along  the  real  f  axis.) 

Second-Order  Processor 


Now  let  x  be  the  sum  of  K  independent  terms  of  the  form  of  (1): 


K 


aks( + 


bklk 


+  cksklk 


+  dksk 


ekV 


k=l 


(8) 


where  real  constants  ak,  bk,  ck,  dk,  ek  can  depend  arbitrarily  on  k, 
and  where  means  msk,  mtk,  standard  deviations  ask,  atk,  and  corre¬ 
lation  coefficients  pk  are  unrestricted  (except  that  ask  _>  0,  a.£k  _>  0, 
|pkl  1  !)•  The  Pair  of  rar|dom  variables  sk,  tk  is  statistically 
independent  of  the  pair  sn,  tn  for  all  k  ^  n.  Thus  random  variable  x  is 
composed  of  a  sum  of  K  groups  of  random  variables,  where  each  group  is 
statistically  independent  of  every  other  group,  but  each  group  itself 
contains  two  statistically  dependent  random  variables. 


This  processor  in  (8)  is  the  general  form  of  interest  in  this  section. 
It  can  be  time-varying  when  the  weights  [ak,  bk,  cR,  dk,  ek]  vary 
with  k,  and  nonstationary  when  the  statistical  parameters  [msk,  mtk, 

CTsk’  atk’  pkl  vary  with  k* 

The  characteristic  function  of  (8)  follows  from  (5)  as 


fx<*>  = 


K 


TT  { i  -  w0  -  ?2°2(k)} 


exp 


k=l 

K 


-1/2 


NQ(k)  -  ffN^(k)  -  f2N2(k) 
3  1  -  i?D1(k)  -  f2D2(k) 


(9) 


*  These  means  can  be  interpreted  as  the  deterministic  signal  components  of 
the  channels  s  and  t,  if  desired. 
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where  the  identification  of  D^k),  etc.,  is  obvious  from  (6).  Only  one 
(continuous)  square  root  and  one  exponential  per  f  value  is  required  in  (9), 
regardless  of  the  number  of  terms  added,  K.  The  mean  and  variance  of  random 
variable  x  in  (8)  follows  from  (9)  as 


(10) 


Any  analytical  attempt  at  determining  the  probability  density  function  or 
cumulative  distribution  function  corresponding  to  characteristic  function  (9) 
would  be  a  formidable  task  indeed.  However,  it  is  a  very  simple  task  via  the 
method  of  [1]  to  get  accurate  numerical  values  for  the  cumulative  and 
exceedance  distribution  functions.  The  program  listing  in  appendix  A 
accomplishes  this  task,  based  upon  characteristic  function  (9)  and  the 
constants  listed  in  (6).  All  the  weights  {ak,  bk,  ck,  dk,  ek]  and 
statistical  parameters  ^msk ,  m^k,  ask,  a^k,  p^are  arbitrary.  Observe 
that  (9)  is  far  more  general  than  the  characteristic  function  considered  in 
[3,  eq.  5],  which  itself  required  a  very  lengthy  analytic  treatment  to  get  the 
probability  density  function  and  cumulative  distribution  function.  In  fact, 
there  is  little  hope  of  getting  any  tractable  analytic  results  for  (9)  when  K 
is  greater  than  2. 

Special  Case  1 

Suppose  weightings  a,  b,  c,  d,  e  in  (8)  are  independent  of  k  and  that 
statistics  as,  a^.,  p  are  also  independent  of  k.  The  decision  variable  x 
in  (8)  then  simplifies  to 


K 


(ID 


k=l 
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Then  D^,  Dg,  ^  are  independent  of  k.  If  we  define  mean  parameters 


K 

M20  =  2  msk  ’  M 
k=l 

K 


K 


02  =  2 
k=l 

K 

M10  =  2  msk  ’  M01  =  2 


mtk  ’  M11 


K 

2 

k=l 


msk  mtk 


"tk 


k=l 


k=l 


(12) 


the  characteristic  function  of  x  in  (9)  then  takes  the  simpler  form 


/  2  VK/2 

yy)  =  fl  -  1?DX  -  fDz)  exP 


Jo  -  igi  -  rN 2 

1  -  i$D,  -  ?2D 


(13) 


where  Dx  and  D2  are  still  given  by  (6),  and 


N0  aM20  +  bM02  +  cMll  +  dM10  +  eM01  ’ 


'1  *  <4ab  -  °tM20  +  1  “sM02  - 


+  (2ae  -  cd)as(a$M01  -  potM1Q)  + 


+  ( 2bd  -  ce)at(otM10  -  pc,.!^)  “ 


v  / 1  .2  2  .  1  2  2.1  . 

-  K  [j  d  0s  +  j  e  at  +  depasat)  , 


N2  =  -  K(ae^  +  bd2  -  cde)(l  -  p2)a2a2 


(14) 


(The  choice  of  K  =  2  and  Ng  =  0  in  (13)  corresponds  to  the  form  given  in 

[3,  eq.  5].)  Observe  that  the  characteristic  function  in  (13)  (and  therefore 
the  performance)  of  the  processor  in  (11)  depends  on  the  means  [my  and 

[my  only  through  the  parameters  [m.  defined  in  (12).  The  mean  and  variance 
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of  random  variable  x  in  (11)  follow  from  characteristic  function  (13)  as 


u 


x 


KDp 


+  2NoDi  - 


KD2  -  2Nj  . 


(15) 


Special  Case  2 


Let  us  also  assume  d  =  0,  e  =  0  in  (11)  above;  then  the  pertinent 
decision  variable  is  given  by 


-2 

k=l 


(ask  +  btk  +  csklk> 


(16) 


and  D2  are  still  given  by  (6),  and  there  follows  from  (14), 


N  =  aM  ■+■  bM  +  rM 
0  20  02  11 


»;.(4ab-c2)(^02M2o.l„2Mo2.pc,s„tHii)  _ 


=  0 


(17) 


The  characteristic  function  of  x  is  given  by  (13),  with  N2  =  0.  The  mean 
and  variance  of  x  in  (16)  are  given  by  (15). 


Fading  for  Special  Case  2 


Let  the  mean  parameters  in  (12)  be  subject  to  slow  fading;  i.e., 

replace 


M20  by  rM20’  M02  by  rM02’  M11  by  rMll’ 


(18) 


where  power  scale  factor  r  has  probability  density  function 
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v  , 

/  \  V  v— i  -vU 

pr(u)  =  ttjt  u  6 


for  u  >  0,  v  >  0; 


r  =  CTr  =  v’  ^(n)  =  for  n  1  (19) 

V 

This  form  of  fading  is  encountered  in  diversity  combination  receivers;  see, 
for  example,  [4,  eq.  9  et  seq.]  and  [5,  eq.  24  et  seq.].  Then  (13),  (17),  and 
(18)  yield  the  conditional  characteristic  function,  for  a  specified  r,  as 


,  f  2  VK/2 

fx(j|r)  =  (1  -  i?D1  exp 


ifr 


-  15N, 


1  -  i£D,  -  f D, 


(20) 


Weighting  (20)  according  to  the  probability  density  function  in  (19),  and 
performing  the  integral,  there  follows,  for  the  characteristic  function  of  the 
decision  variable  x  in  (16),  the  result 


_  T1  -  w i  -  ?2°2Y " 2 

1  -  1?(DX  +  Nq/v)  -f2(D2 


N,/v 


(21) 


(The  limit  of  (21)  as  v-*+»is  again  (13)  with  =  0,  as  in  (17);  this 
agrees  with  the  fact  that  the  corresponding  limit  of  the  probability  density 
function  in  (19)  is  Pr(u)  =  S(u-l).)  The  mean  and  variance  of  x  in  (16) 
follow  from  characteristic  function  (21)  as 


=  \  KD2  +  2NqD1  -  KD2  -  2N^  +  Nq2/ v  .  (22) 

Observe  that  mean  ux  is  independent  of  v,  the  power  law  in  fading  (19).  A 
program  for  the  cumulative  and  exceedance  distribution  functions  corresponding 
to  characteristic  function  (21)  is  given  in  appendix  B. 
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SPECIAL  FORMS  OF  SECOND-ORDER  PROCESSOR  (8) 


Before  embarking  on  the  analysis  of  the  other  two  classes  of  processors, 
we  will  explicitly  detail  some  of  the  special  forms  that  processor  (8)  reduces 
to,  under  particular  selections  of  the  weightings  and  statistical  parameters. 

A  rather  broad  collection  of  typical  processors  will  be  seen  to  be  included. 

In  the  following,  any  unspecified  weights  £ak,  bk,  ck,  dk,  ek]  are  zero, 
and  any  unspecified  statistical  parameters  that  do  not  appear  in  the  final 
characteristic  function  are  irrelevant. 

I.  Gaussian 


K 


k=l 


K  K 


V?>  -  exe  msk  -kf  2 


2 


k=l  k=l 


II.  Chi-square  of  K  Degrees  of  Freedom 


ak  *  l'  msk  *  °>  ”sk  *  1 


K 


k=l 


fx(?)  =  (l  -  i?2)  K/2 
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III.  Non-Central  Chi-Square  (QM  Distribution  if  K  =  2M) 

Tsk  =  as 


\  =  1>  oc>.  =  a. 


K 


k=l 


fx(f)  =  (l  -  1J2 »0 


— K  /  2 


exp 


K 

if^  m 

k=l 


sk 


1  -  1f2a‘ 


IV.  Weighted  Energy  Detector 


ak  ^  0,  dk  =£  0 


■I 


(aksk  +  dksk} 


k.l 


U?>  = 


j  {i  -  *2v4 

k=l 


1 

2 


exp 


akmsk  +  dkmsk  +  ^Kask 


K  SK  K  SK 

k=i  1  '  *2V 


k  sk 


V.  Weighted  Cross-Correl ator 


ck  ^  0 


K 

x '  2ckVk 
k=l 

D^k)  =  2ckpkaskatk, 

D2(k)  =  -C2(l  -  Pk)askatk, 
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N0(k)  "  Ckmskmtk’ 

Nl(k)  =  -ck(-2  msk  atk  +  2  mtkask  "  pkmskmtkCTskCTtk^ * 


N2(k) 


=  0 


Characteristic  function  f  (£)  is  given  by  (9). 


VI.  Two-Channel  Energy  Detector 

ak  ^  0,  bk  ^  0 

K 

x  =  2  (aksk  +  bktk^ 
k=l 


D1(k) 

02(k) 

N0(k) 

Nx(k) 

N2(k) 


2(Vsk  +  bkCTtk}  » 


4akbk^  “  » 


k'  sk  tk 


akmsk  +  bkmtk’ 


4akV?  mskCTtk  +  2  mtkask  “  pkmskmtkCTskatk^’ 


0. 


Characteristic  function  fx(£)  is  given  by  (9).  A  simple  application  of  this 
particular  processor  was  encountered  in  [6,  eqs.  25-26]. 
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VII.  Two-Channel  Energy  Detector  and  Cross-Correlator 


ak  ^  0,  4=  0,  C|<  ±  0 


=  1 

k=l 


(aksk  +  bkbk  *  ckskV 


Dl(k>  ‘  2(Vsk  *  Vtk  +  'kV’sk'tk* 

02(k)  -  (4akbk  -  c2)(l  -  p2),^2 

N0(k)  ■  Vsk  *  bkratk  +  ckmskmtk  • 


Nl(k)  -  (4akbk  -  Ck)(I  ™sk<’tk  +  \  "tk'sk  -  Skmsk”tk''5k,’tk> 


M2(k) 


=  0 


Characteristic  function  f  (F)  i s  given  by  (9) . 
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NARROWBAND  CROSS-CORRELATOR 


s^t) 


s9(t) 


+  n^t) 


n2(t) 


\ _ * 

Lowpass 

z(t)  > 

Sample,  Weight,  . 

V  ^ 

Filter 

and  Accumulate 

— -  > 

Figure  1.  Narrowband  Cross-Correl ator 


The  processor  of  interest  in  this  section  is  depicted  in  figure  1.  Input 
signals  s^(t)  and  s2(t)  are  arbitrary  deterministic  narrowband  real  waveforms: 

s j ( t )  =  R&[sj(t)  exp(  i2irfQt)J  =  A ^ ( t )  cos(2TrfQt  +  Pj(t))  = 

=  aj(t)  cos ( 2irf Qt )  -  bj ( t )  sin(2irf0t)  for  j  =  1,2,  (23) 

where  input  signal  complex  envelope 


Sj(t)  =  A j ( t )  exp( iPj(t) )  =  a j ( t )  +  ibj(t)  (24) 

in  terms  of  polar  or  rectangular  low-frequency  components,  respectively. 

Input  noises  n^t)  and  n2(t)  are  zero-mean  correlated  narrowband 
jointly-Gaussian  processes  which  may  be  nonstationary: 

n j ( t )  =  Re[n^(t)  exp(i2rfQt)]  =  x^(t)  cos(2rfQt)  -  y.{t)  sin(2irfQt),  (25) 
where  noise  complex  envelope 

J2j(t)  =  Xj(t)  +  i y j ( t )  for  j  =  1,2.  (26) 

The  statistics  of  the  input  noise  complex  envelopes  are  arbitrary: 
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n~Tt]  jt^C t )  =  0  for  all  j,  in. 


(27) 


The  quantities  a2,  y  can  all  vary  with  time  t,  for  nonstationary  noise 
processes.  There  follows,  for  the  statistics  of  the  in-phase  and  quadrature 
components  defined  in  (25), 


2  2  2  - 

X1  =  yl  =  ffl*  xlyl  =  0  ’ 


X1X2  '  yly2  "  ct1ct2p  » 


(28) 


x2yi  =  =  ^1^2X  • 


The  reason  for  breaking  out  this  narrowband  cross-correl ator  as  a 
separate  problem  is  now  apparent  from  (28).  Namely,  at  each  time  instant,  a 
group  of  four  random  variables  are  statistically  dependent  on  each  other. 

This  case  does  not  fall  into  the  framework  of  (8)  above,  since  only  two  random 
variables  were  dependent  there. 

Using  the  narrowband  character  of  all  the  waveforms  in  (24)  and  (26),  the 
lowpass  filter  output  in  figure  1  may  be  expressed  as 


z(t)  =  -|[x1(t)  +  a1(t)][x2(t)  +  a2 ( t ) ]  +|[y1(t)  +  b1(t)][y2(t)  +  b2(t)].  (29) 
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The  final  system  output  in  figure  1  is  the  weighted  sum  of  K  terms, 

K 

v-5.  w(k)  z(tk)  ,  (30) 

k=l 


where  it  is  assumed  that  the  time  separations  between  samples  at  instants 
ft^  lead  to  statistically  independent  random  variables  ^(t^)}  .  The 
weights  and  statistics  can  change  with  sample  time  t^,  in  an  arbitrary 
fashion. 

Based  upon  the  method  in  [7],  we  find  the  characteristic  function  of  z(t) 
in  (29)  to  be  given  by 


fz(T.t)  = 


i-i?o1  +  r d2 


exp 


1?  "o  +  (5N1 


1-ifDj 


where  the  constants  (in  their  most  compact  form)  are  given  by 


D1  =  ai°2p  ’ 


n  _  1  2  2n  2 


(31) 


M0  _  2 ^ ala2  +  blb2^  ’ 


N 


1  "  8l 2^al  +  bl^  +  CTl^a2  +  b2^  "  201a2p(a1a2  +  b^)  -  20^(3^  -  a^^J . 


(The  characteristic  function  and  constants  in  (31)  and  (32)  are  not  to  be 
interchanged  or  confused  with  any  earlier  results  in  previous  sections.  In 
fact,  observe  there  is  no  square  root  involved  in  (31).)  All  of  the 
parameters  in  (32)  can  vary  with  time  t. 


(32) 


In  terms  of  the  signal  polar  definitions  in  (24)  and  the  complex  noise 
correlation  coefficient  y  in  (27),  alternative  expressions  to  (32)  (where  we 
have  emphasized  the  t-dependence)  are 
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Dj_  =  cr1  ( t )  o2( t)  Re[y(t)}  =  ffx(t)  o2(t)  | Y ( t )J  COS  <i){t)  , 


D2  =  \  a2(t)  o\{ t)  (l-|Y(t)|  2)  , 

Nq  =  \  Re|^(t)  s2(t)}  =  J  A^(t)  A2(t)  cos[P1(t)-P2(t)]  , 

N1  =  ^£2^  2  +  al^J-2^  2-2  al^t^  Re{li(t)  I^)  y(t)jj  = 


=  A  a2(t)  Aj(t)  +  Oj(t)  A2(t)  -  2  a^t)  o2(t)  A^t)  A2(t)  |y(t)j  COsC P ( t ) -P 2( t ) -(6 ( t )  1 1 
L  (33)  J 

The  mean  and  variance  of  z(t)  in  (29)  follow  from  (31)  as 


uz  =  D1  +  N0  ’ 


“z  =  Dl  +  2“2  +  20iN0  *  2N1 


(34) 


Finally,  the  characteristic  function  of  the  narrowband  cross-correl ator 
output  v  in  (30)  follows  from  (31)  as 


fv<F>  -  IT  fz(?“<k)'tk)  = 
k=l 


"fj“  {l-ilw(k)  Dx(k)  +f2w2(k)  D2(k)l 
k=l  1 


-1 


exp 


K  w(k)  .Uk)  +ijw2(k)  Nl(k)n 
i  |  \  - r _ ± _ 

j“i  l-i>(k)  Dx(k)  +  |2w2(k)  D2(k) 


where  we  have  allowed  all  the  parameters  in  (32)  and  (33)  to  vary  with  time 
tk.  The  mean  and  variance  of  output  v  follow  from  (35)  as 


(35) 


uv  = 


w(k)  [D^k)  +  NQ(k)], 


k=l 


%  =  w2(k)  fDl(k)  +  2D2(k)  +  2DiO<)N0(k)  +  2N1(k)7.  (36) 

k=l  J 
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A  program  for  the  evaluation  of  the  cumulative  and  exceedance  distribution 
functions  via  (35)  is  given  in  appendix  C. 

In  comparison  with  earlier  results  in  [7],  we  have  obtained  the  following 
extensions  here: 

1.  The  input  signals  are  arbitrary  narrowband  waveforms;  they  are  not 
limited  to  two  sine  waves  at  the  same  frequency; 

2.  The  Gaussian  input  noises  can  be  nonstationary; 

3.  The  number  of  terms  summed  to  yield  the  narrowband  cross-correlator 
output  can  be  greater  than  1; 

4.  The  characteristic  function  is  in  its  most  compact  form,  and  the 
constants  are  expressed  directly  in  terms  of  given  quantities, 
having  eliminated  all  auxiliary  variables. 

Output  Signal-to-Noise  Ratio 

It  is  sometimes  desirable  to  have  simple  expressions  for  the  output 
signal-to-noise  ratio  of  the  narrowband  cross-correl ator  in  figure  1.  In 
terms  of  the  lowpass  filter  output  z(t),  we  observe  first  from  ( 32) — ( 34)  that 

uz(s)  =  uz(s+n)  -  uz(n)  =  Nq  =  j  A1  A2  cos^-P^  .  (37) 

We  then  have  two  alternative  definitions  of  the  signal-to-noise  ratio  at  the 
lowpass  filter  output: 


Rz(n) 


uz(s)  A^  A^  cos2(P^-P2) 

"TTT  =  “  2  2 777  2  2s  ’ 

CTz(n)  2  cr 2 ( 1+P  -x  ) 


Rz(s+n) 


uz(s)  A2  A2  cos2(PrP2) 

a2(s+n)  4(D2  +  2D2  +  2D1N0  +  2N^) 


(38) 


18 


TR  7035 


These  closed  form  expressions  allow  for  arbitrary  noise  correlations  and  are 
considerably  simpler  than  [7,  eqs.  41-43].  The  signal-to-noise  ratios  of 
system  output  v  in  figure  1  are  K  times  greater  than  either  form  in  (38). 

Specialization  to  Narrowband  Energy  Detector 

If  the  signal  and  noise  parameters  in  (24)  and  (27)  are  chosen  as 

ax(t)  =  a2(t)  -  a(t)  , 

bx(t)  -  b2(t)  =  b(t)  , 

a^(t)  =  a2(t)  =  a(t)  , 

p(t)  =  1,  X(t)  =  0  ,  (39) 

then  figure  1  reduces  to  identical  input  channels,  that  is,  a  narrowband 
energy  detector.  There  follows  from  (32), 


D,  =  cr  ( t ) ,  D9  =  0  , 


«  ^a2(t)  +  b2(t))  = 


\  A2(t),  Nj_  =  0  , 


(40) 


and  (31)  becomes 


f7(?,t)  =  - i-. - 

z  i  -  iVu) 


exp 


IS- 


A2(t)/2 


(41) 


1  -  if*  (t) 

Corresponding  results  for  the  system  output  v  are  easily  obtained  from  this. 


19 


TR  7035 


REDUCTION  OF  HERMITIAN  AND  LINEAR  FORM 

The  most  general  case  of  interest  in  this  section  is  as  follows: 
random  complex  matrix 


X  =  [x1  x2  ...  xM] 


(42) 


is  Mxl;  constant  complex  matrix 


A  -  [a1  a2  ...  aM] 


(43) 


is  Mxl;  and  constant  complex  matrix 


B  = 


bll  '*•  blM 


bMl  *'•  bh 


(44) 


is  MxM  and  Hermitian.  The  Hermitian  and  linear  form  we  consider  is 


q  =  XHBX  +  -|(XHA+AHX)  = 

M  M 

*  5  xmbmnxn  +  7  X  (xmVVm>'  <45> 

m,n=l  m=l 

which  is  real.  Random  variable  q  is  a  weighted  sum  of  all  possible  products 

°f  fxm"j  and  [xJ.  plus  linear  combinations.*  A  and  B  are  called  the 
weighting  matrices. 


*  For  M=2  or  4,  and  real  variables  and  weights,  (45)  reduces  to  the  earlier 
forms  given  in  (1)  and  (29). 
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We  will  concentrate  in  this  section  on  reducing  form  (45)  to  a  weighted 
sum  of  squares  of  uncorrelated  random  variables.  This  stepping  stone  does  not 
require  any  Gaussian  assumptions  on  X  and  is  therefore  useful  as  a  separate 
item. 


The  relevant  statistics  pertaining  to  random  vector  X  are 


X  =  E  (mean  matrix) , 
X  =  X-  X  =  X-E, 


Cov[x)  =  XXH  =  K  (covariance  matrix),  (46) 

where  statistics  matrices  E  and  K  are  given.  MxM  matrix  K  is  always  Hermitian 
and  non-negative  definite.  We  assume  K  is  positive  definite;  otherwise 
eliminate  the  linearly  dependent  components  of  X.  We  allow  xm  and  xn  to 
be  correlated  with  each  other  for  any  m  and  n;  this  situation  is  much  more 
general  than  the  investigations  above. 

Let  C  be  a  constant  MxM  matrix  and  form  the  linearly  transformed  variables 

W  =  CHX  =  [w1  w2  ...  wm]T.  (47) 

Then  the  statistics  of  W  are  given  by 

W  =  CHE, 

W  =  w  -  w  =  c"x. 


Covfw}  =  wwH  =  CHX?C  =  CHkc. 
Also,  from  (47),  since 


(48) 


X  =  C"HW, 


(49) 
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then  we  can  express  (45)  as 

q  =  WHC_1BC"HW  +  -j(WHD  +  DHW), 

where  we  define  constant  Mxl  matrix 

D  =  C-1A  =  [d1  d2  ...  dM]T. 

We  want  to  have,  from  (48)  and  (50), 

CHKC  =  I 


(50) 


(51) 


(52) 


and 


C  1BC  H  =_A.=  diag(x1  x2  ...  xM) ;  i.e.  CHB  1C  =_/L  1  , 
for  then,  in  addition  to  the  relation  between  the  means, 

¥  =  CHE, 


(53) 


(54) 


we  have  the  desirable  properties 


Cov{w]  =  I, 


(55) 


and 


U  1  U  LI  M 

q  =  Wn4-W  +  j(WnD  +  DnW)  =  ^  xmlwm 

m=l  m  m 


+  Re 


m=l 


* 

dmwm 
m  m 


(56) 


That  is,  the  random  vector  W  given  by  (47)  is  composed  of  uncorrelated  unit- 
variance  components,  and  q  is  a  weighted  sum  of  magnitude-squares  of  these 
components,  in  addition  to  a  linear  sum. 


We  now  have  to  address  the  problem  of  determining  the  MxM  matrices  C  and 
_A.in  (52)  and  (53).  From  [8,  p.  106,  Theorem  2],  we  identify 
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M— K,  K-*  B_1,  (57) 

then  according  to  [8,  p.  107,  eq.  29],  we  must  solve  for  C  and  A  in  the 
equation 


B-1C  =  KCA-1,  i.e.  BKC  =  CA.  (58) 

So  the  only  matrix  that  need  be  considered  is  the  MxM  product  BK.  C  is  the 
modal  matrix,  and  A  the  eigenvalue  matrix,  of  BK.  Also,  from  (51), 

D  =  CHA,  since  C_1  =  CH.  (59) 

Letting  C  =  [C^  ...  C^],  where  eigenvector  is  a  Mxl 

matrix,  (58)  can  be  expressed  as 

BKC^  =  xmC(m)  for  1  <_  m  <_  M.  (60) 

Several  important  properties  hold  for  A  and  C: 

The  are  all  real,  but  can  be  positive,  zero,  or  negative. 

If  K  and  B  are  real,  then  C  is  real.  (61) 

If  B  is  positive  definite,  then  xm  >  0  for  1  <  m  <  M. 

If  A  =  0  and  E  =  0,  there  is  no  need  to  solve  (58)  for  C, 
because  D  =  0  and  W  =  0. 
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QUADRATIC  AND  LINEAR  FORM 

If  random  vector  X  is  real  Gaussian,  if  A  is  real,  and  if  B  is  real 
symmetric,  then  mean  E  and  covariance  K  are  real,  and  it  follows  that  modal 
matrix  C  is  also  real.  Also  from  (47)  and  (59),  W  and  D  are  real.  Equation 
(45)  reduces  to 

XT  JL  M 

q  =  X'BX  +  x'a  =  \  xb  x  +  ^\ax  ,  (62) 

^  m  mn  n  ^  mm’  v  / 

m,n=l  m=l 

which  is  a  quadratic  form  and  linear  form. 


Letting  mean  W  in  (54)  be  expressed  as 


W  -  [v1  \>2  •  •  •  > 


(63) 


the  Gaussian  character  of  X  and  the  linear  transformation  (47)  allow  us  to 
write  the  probability  density  function  of  W  as  a  product: 


P(W) 


M  _ 

=  II  f(2ir)"1/2  6XP(l  (wm  "  vm)2)] 


(64) 


Here  we  used  property  (55).  Since  we  now  have,  from  (56), 


M 


q  “  f-  <Vm  +  drawra>- 

m=l 


(65) 


the  characteristic  function  of  q  is 


fq(?)  =  exp(i^q)  =  exp ( i 


+  d",“ra)) 


n 

TT 

m=l 


-1/2  r 


exp 


'52 


M  x  v2  +  d  v  +  i£d2/2 
mm  mm  >  m 


m=l 


l2xm^ 


,  (66) 
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where  the  square  root  must  be  a  continuous  function  off,  not  a  principal 
value  square  root.*  Notice  that  only  one  square  root  and  one  exponential  is 
required  per  f  value.  Observe  that  the  characteristic  function  depends  on  the 

separate  values  ^  and  fd^p  not  merely  on  their  sums.  If  A  =  E  =  0,  the 

exponential  is  unity,  by  virtue  of  (54)  and  (59).  And  if  M=2,  (66)  reduces  to 

(5),  while  M=4  leads  to  form  (31). 

To  summarize,  the  characteristic  function  f  (f)  in  (66)  for  random 
variable  q  in  (62)  requires  the  constants  [xj ,  [dj,  and  [vm]  for  1  <  m  <  M. 
The  initially  given  quantities  are  weighting  matrices  A,  B  and  statistics 
matrices  E,  K.  We  first  solve  the  equation  (58), 

BKC  =  CA,  (67) 

for  eigenvalue  matrix  JL  and  modal  matrix  C  corresponding  to  BK.  Then 

A  =  diag(x:  x2  ...  xM)  , 

D  =  CTA  =  [d:  d2  ...  dM]T  , 

W  =  CTE  =  [Vl  \>2  ...  vm]T  .  (68) 

If  the  mean  of  input  X  is  zero,  E  =  0,  and  if  the  linear  weighting  form  is 
zero,  A  =  0,  then  there  is  no  need  to  solve  for  modal  matrix  C  of  BK  in 

(67).  Then  D  =  W  =  0  and  the  exponential  term  in  (66)  is  unity.  One  only 
need  compute  eigenvalue  matrixAof  BK  in  this  case. 

A  program  for  the  evaluation  of  the  cumulative  and  exceedance 
distribution  functions  corresponding  to  characteristic  function  (66)  is  listed 
in  appendix  D.  The  inputs  to  the  program  are  considered  to  be  M,  fxA, 


*  That  is,  the  square  root  is  the  analytic  continuation  of  the  function 
defined  as  1  at  f  =0. 
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[dm}>  [vm];  that  is,  it  is  presumed  that  (67)  and  (68)  have  already  been 
solved  prior  to  use  of  the  program. 

The  cumulants  of  q  are  obtained  from  (66)  as 


*q(n)  = 


^  (x  +  X  v  +  d  v  )  =  u 
™  rn  m  mm  q 


for  n  =  1 


i 


M  n  9 

9n-l,  1U  ^  xn  c 
2  ( n— 1)1  2L  m 


V_ 


m=l 


x2  +  n(x  v  +id') 
m  \  m  m  2  ny 


for  n  >  2 


(69) 


In  particular,  the  variance  of  q  is 


M 

Xa<2>  -  2^ 

q  m=l 


+  2 


(70) 


If  another  random  variable  is  formed  by  the  sum  of  several  independent 

random  variables  q.  with  the  form  (62),  but  with  different  sizes  M,,  the 
J  J 

new  characteristic  function  is  the  product  of  terms  like  (66). 


Breakdown  of  X  into  Two  Components 


It  is  useful  to  investigate  a  particular  version  of  the  general  results 
above,  because  the  resultant  forms  correspond  to  some  often-realized  practical 
energy  detectors  and  correlators.  We  let  M  =  2N,  and 


u 

'b11  b12~ 

(a'1*! 

II 

X 

V 

_  -i 

,  B  = 

J 

B21  B22 

II 

<c 

I - 

no 

where  U ,  V,  are  Nxl  real  matrices,  and  are  NxN  real 

matrices.  Also  and  B22  are  symmetric,  while  B2-'-  =  B-*-2  .  Then 
(62)  can  be  expressed  as 
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Hi1  vt1 

V1  B12" 

u 

+  r„T  vTl 

r^i 

Lu  VJ 

1 

_821  B22. 

; 

V 

-  - 

+  l_U  V  J 

CM 

_ 1 

q  =  XTBX  +  XTA 


=  UTBnU  +  UTB12V  +  VTB21U  +  VTB22V  +  uV1)  +  vV2>  = 

=  utbuu  +  2utb12v  +  vtb22v  +  uV1)  +  vta(2>  = 


N 


XT  /  ,11  .  .  ,12 

=  ' >  (u„  b +  2u„  b v 


m  mn  n  m  mn 


+  v  b22  v  \  +  fu 
n  m  mn  n)  n  n 

7  n=l  x 


(D  +  v 

n  n  j 


m,n=l  N 

=  all  possible  auto  and  cross  combinations  of  random 
variables  ju^  and  fv^jj*  ,  plus  linear  combinations. 

We  also  have,  from  (47)  and  (68), 


W  =  CTX  =  CT 


W  =  C 


=  C1 


D  =  C1 


,(2) 


QTl 

|f 

|t  jf] 

PJ 

LL 

1  M 

K  =  Cov{x]  =  XXT  = 

Then  the  fundamental  matrix  required  in  (67)  is  expressable  as 


K  K 

uu  uv 

K  K 

vu  vv 


V1 

B12 

1 

\u 

K 

uv 

B21 

b22_ 

‘  ; 

K 

_  vu 

- 1 

> 

> 

BK  = 


which  is  a  2Nx2N  matrix.  Also  random  variable  (65)  is  now 


2N 

q  =  ''S  (aw)"  +  dw 
\  mm  mm 

m=l 


(72) 


(73) 


(74) 


(75) 
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which  has  2N  terms.  The  characteristic  function  of  q,  in  general,  follows 
from  (66)  and  (68)  as 


2N 

-1/2 

2N 

O  O 

IT  f1-121/} 

exp 

x  +  d  v  +  iTd„/2 

mm  mm  >  m 

1-iZx  f 

jn=l  _ 

, 

m=l 

III  9 

(76) 


where 


W  =  [v1  . . . 


V2N^ 


D  =  [dj_  d2 


"A(lf 

a(2) 


(77) 


(If  A^'  =  A^  =  Eu  =  Ev  =  0,  then  D  =  0  and  W  =  0,  and  there  is  no 

need  to  solve  for  modal  matrix  C;  the  exponential  in  f  (£■)  in  (76)  is  then 
unity. ) 


As  a  special  case,  if  A  =  0,  B11  =  0,  B22  =  0,  then  (71)  and  (73) 

yield 


and  (72)  gives 


B  = 


0  B12' 

B21  0 


D  =  0  , 


T  12 

q  =  2U  B1  v  =  2 


m,n=l 


.12 

um  bm„  v„  = 
m  mn  n 


(78) 


=  all  possible  cross  combinations  of 


Then  (74)  specializes  to 


0 

B12" 

K 

K 

b12k 

b12k 

BK  = 

uu 

uv 

— 

VU 

vv 

B21 

0 

1 

K 

K 

b21k 

b21k 

— 

# 

L  vu 

vvj 

L  uu 

UV 

(79) 


(80) 
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and  (75)  reduces  to 


2N 


q  = 


x_w2 
m  m 


m=l 


with  characteristic  function 


~2N  r  -> 

-1/2 

2N  XmVm 

V?)  - 

ir 1 1  -  i2i4 

exp 

I 

m=l 

L  m=1  J 

following  directly  from  (76)  and  (78). 


(81) 


(82) 


For  the  particular  example  of 


B12  =  J  diagO^  £2 - GN)  +  \  [g-j.  92  •••  9N]T  [h1  h2  ...  hN],  (83) 


then 


N  \  /  N 


"  '  ^Tl  "i"U"Vn  +  S  9nU"J  (  n?i  h"Vn  1  ’ 


(84) 


with  the  same  characteristic  function  (82). 


12  1 

As  a  still  more-special  case,  let  B  =  I;  then  (79)  and  (81)  give  the 


simple  cross-correlator  (but  with  correlated  inputs  for  all  time  separations) 

__  .  (85) 

=1  m=l 

and  (80)  and  (77)  become 


N  2N 

“  -  j>  v„  - 

n=l  m=l 


x  w2 
m  m 


BK  =  j 


K 

K 

Fe  ~1 

vu 

vv 

s=l 

ii 

o 

— H 

u 

K 

K 

E 

uu 

uv_ 

v_ 

=  c^ 


1 


V2N^ 


(86) 
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The  important  equations  that  must  be  solved  are  always 

BKC  =  CJL  (87) 


or 

BKC^  =  xmC(m)  for  1  £  m  £  M  =  2N,  (88) 

where  all  matrices  B,  K,  C,j\_  are  2Nx2N.  The  characteristic  function  of 
(85)  is  again  (82). 

Special  Case  of  Correlator  (85) 


Here  we  let  components  U  and  V  have  the  same  covariance  and  a  scaled 
cross-correl ation ;  that  is,  let 


K 

uu 


K  =  K  , 

VV  0 


K  =  PK 

UV  0 


(89) 


where  p  is  a  scale  factor.  This  case  corresponds,  for  example,  to  a  common 
signal  in  two  independent  components:' 


u(t)  =  s ( t )  +  n1(t), 
v ( t )  =  s ( t )  +  n2(t). 


(90) 


where  s(t),  n^(t),  n£(t)  are  all  independent  and  have  a  common  covariance. 
Then  (86)  becomes 


BK 


(91) 


Now  suppose  that  we  can  determine  the  NxN  eigenvalue  matrix  f  and  modal 
matrix  Q  of  Kqj  that  is 
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K0Q  =  or  ;  r=  diag(Yl>  y2  •••  yn). 
Then  we  have  the  standard  relations  [8] 

K0  =  QTQT  where  QQT  =  I. 

We  can  now  express  the  2Nx2N  matrix  in  (91)  as 


BK  =  4 


QPrQT 

orqT 

_qpqt 

QproTj 

Q  o 

0  Q 

and  2Nx2N  identity  matrix 

fQ  cPr 


JP? 

\p? 


qt  0 


0  Q1 


I 


2N 


0  Q 


I  0 

0  I 


'  T 

Q1  0 


0  Q1 


There  follows 


BK  -  x  1 - 

But  the  middle  matrix  in 
form 


Q 

o1 

|»r  -  xi 

?r  " 

"  T  ~1 

Q1  0 

0 

Q 

0  qt 
■.  — < 

(96)  can  be  developed  in  detail  in  the 


(92) 

(93) 

(94) 

(95) 

(96) 

partitioned 
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x 

1 

1 

l 

Vl 

1 

7pyN  "  x| 

1 

Fn 

1 

Vl 

1 

1 

1 

1 

'jPYj.  ~  X 

» 

(97) 

1 

1  ' 
Fn 

i 

Ipyn  ■  x 

This  matrix  is  singular  when  the  kth  row  is  equal  to,  or  the  negative  of,  the 
k+N  tji  row.  This  leads  to  the  eigenvalues  M  ^  of  matrix  BK: 


^1  —  2^p+^y1’  ***’  —  2^p+^yN  ’ 

x[\|+l  =  2'(p-1)y^j  •••»  x2N  =  '  ( 98) 


Thus  we  need  only  solve  for  the  N  eigenvalues  of  matrik  KQ ,  and 

then  use  them  as  above  to  determine  all  2N  eigenvalues  of  BK;  this  is  a 
significant  shortcut. 


If  also  Eu  =  Ey  =  0,  then  W  =  0  from  (86),  and  the  characteristic 
function  of  q  in  (85)  follows  from  (82)  and  (98)  as 


y?>  - 


m=l 


-1/2 


IN 

IT  ?2} 

m=l 


-1/2 


(99) 


This  is  a  generalization  of  [1,  eq.  54],  which  held  for  a  single  pair  of 
Gaussian  random  variables. 
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EXAMPLES 

The  program  listed  in  appendix  A  for  the  second-order  processor  (8)  and 
attendant  characteristic  function  (9)  has  been  employed  to  yield  the  result  in 
figure  2.  The  particular  values  for  the  number  of  terms  K,  the  weights,  and 
the  input  statistics  are  listed  in  lines  20-120.  There  is  no  physical 
significance  attached  to  this  particular  example;  rather  it  has  been  run 
simply  to  illustrate  the  extreme  generality  that  the  technique  is  capable  of. 
Some  negative  values  for  the  weights,  means,  and  correlation  coefficients  have 
been  employed  to  emphasize  this  generality.  This  simple  example  (and  others 
to  follow)  can  be  used  as  a  check  case  on  any  user-written  program  to  evaluate 
cumulative  and  exceedance  distribution  functions. 

The  selection  of  parameters  L,  a,  b  in  lines  130-150  is  discussed  in 
detail  in  [1];  the  reader  is  referred  there  for  the  deleterious  effects  that 
can  occur  for  improper  choices  of  L,  a,  b.  The  selection  of  Mf,  the  FFT 
size  in  line  160,  is  rather  arbitrary;  it  controls  the  spacing  at  which  the 
probability  distributions  are  computed,  but  has  no  effect  upon  the  accuracy  of 
the  results  (except  for  round-off  noise).  Additional  computational  details  on 
the  particular  program  for  characteristic  function  (9)  are  given  in  appendix  A. 

The  ordinate  scale  for  figure  2  is  a  logarithmic  one.  The  lower  right 
end  of  the  exceedance  distribution  function  curve  decreases  smoothly  to  the 
region  IE— 11 ,  where  roundoff  noise  is  encountered.  The  exceedance 
distribution  function  values  continue  to  decrease  with  x  until,  finally, 
negative  values  (due  to  roundoff  noise)  are  generated.  For  negative 
probability  values,  the  logarithm  of  the  absolute  value  is  plotted,  but 
mirrored  below  the  IE— 1 2  level.  These  values  have  no  physical  significance, 
of  course;  they  are  plotted  to  illustrate  the  level  of  accuracy  attainable  by 
this  procedure  with  appropriate  choices  of  L  and  a. 

The  rates  of  decay  of  the  cumulative  and  exceedance  distribution 
functions  in  figure  2  are  markedly  different  for  this  particular  example. 
Additionally,  since  the  decays  are  both  linear  on  this  logarithmic  ordinate, 
it  means  that  both  tail  distributions  are  exponential,  not  Gaussian.  These 
attributes  of  the  cumulative  and  exceedance  distribution  functions  are  easily 
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Figure  2.  Distributions  for  Second-Order  Processor 
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and  quickly  discernible  by  use  of  the  numerical  technique  in  [1],  for  a 
limitless  variety  of  weights  and  input  statistics,  with  a  minimum  of  effort  on 
the  part  of  the  user. 

As  a  check  on  the  program  in  appendix  A,  the  second-order  processor  in 
(8)  was  simulated,  and  10,000  independent  trials  were  used  to  determine  its 
performance  for  the  exact  same  parameters  as  used  for  figure  2  above.  The 
program  is  listed  in  appendix  E  and  the  results  are  given  in  figure  3.  The 
corroboration  is  excellent,  even  near  the  IE— 4  probability  level. 

As  the  number  of  terms,  K,  in  the  second-order  processor  (8)  is 
increased,  and  if  the  statistics  are  identical,  the  random  variable  x  should 
approach  Gaussian,  at  least  near  its  mean.  The  example  in  figure  4  was  run 
for  K  =  10,  and  all  weights  and  statistics  independent  of  k;  the  particular 
choices  were 

a  =  .6,  b  =  -.6,  c  =  .3,  d  =  -.2,  e  =  .2, 

ms  =  -5,  mt  =  -.5,  0$  =  1,  ot  =  1,  p  =  .4,  (100) 

L  -  4,  A  =  .05,  b  =  20ir,  Mf  =  256. 

The  cumulative  and  exceedance  distribution  functions  in  figure  4  both  display 
a  parabolic  shape  near  the  mean  of  x,  which  signifies  Gaussian  behavior  of  the 
random  variable,  as  expected.  However,  on  the  tails,  the  distributions  are 
tending  to  linear,  which  means  an  exponential  decay  there.  This  observation 
for  this  example  confirms  the  comments  of  [7,  p.  673]. 

The  cumulative  and  exceedance  distribution  functions  for  an  example  of 
the  second-order  processor  with  fading  are  displayed  in  figure  5,  as 
determined  from  characteristic  function  (21)  and  the  corresponding  program  in 
appendix  B.  The  power  law,  v,  for  the  fading  probability  density  function 
(19)  is  2.7  for  this  example,  but  can  be  easily  changed.  The  particular 
constants  employed  are  listed  in  lines  20-110  in  appendix  B. 

An  example  of  the  distributions  for  the  narrowband  cross-correl ator  of 
figure  1  is  presented  in  figure  6,  as  evaluated  from  characteristic  function 
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-2&7T  -/Ott  0  x  IQrr  20ir 


Figure  4.  Second-Order  Processor,  K=10,  Identical  Statistics 
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Figure  5.  Second-Order  Processor  with  Fading 
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Figure  B.  Distributions  for  Narrouband  Cross-Cor re  1 ator 
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(31)  and  the  program  in  appendix  C.  The  weightings,  signal  components,  and 
noise  statistics  have  no  special  values  or  interrelationships;  the  particular 
values  used  here  are  listed  in  lines  20-110. 

The  distributions  for  the  reduced  quadratic  and  linear  form  (65)  and 
accompanying  characteristic  function  (66)  are  presented  in  figure  7  for  the 
numerical  example  employed  in  the  program  listing  in  appendix  D.  If  the  given 
form  is  instead  that  of  (62),  then  ( 67)— ( 68)  must  first  be  solved  before  the 
program  in  appendix  D  can  be  employed;  that  is,  one  must  augment  these  results 
with  the  capability  for  extracting  the  eigenvalues  (and  eigenvectors  in  some 
cases)  of  the  MxM  matrix  BK.  The  size  of  the  FFT,  has  been  increased  to 
1024  in  figure  7;  this  results  in  finer  spacing  of  the  distribution  values  and 
additional  spikes  in  the  round-off  noise  region  centered  about  IE— 12 . 
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Figure  7.  Distributions  -for  Quadratic  and  Linear  Form 


41 


TR  7035 


SUMMARY  AND  DISCUSSION 


Closed  form  expressions  for  the  characteristic  functions  of  the  decision 
variables  of  three  classes  of  second-order  processors  have  been  derived.  The 
input  noise  to  the  processors  must  be  Gaussian,  but  it  can  be  nonstationary 
with  arbitrary  statistics.  Programs  for  the  direct  evaluation  of  the  exact 
cumulative  and  exceedance  distribution  functions  have  been  generated  and  then 
exercised  for  completely  general  values  of  the  weights,  signal  parameters,  and 
noise  statistics.  There  is  no  assumption  needed  about  a  large  number  of 
statistically  independent  contributors,  nor  need  any  signal-to-noise  ratio  be 
either  small  or  large.  The  first  two  classes  of  processors  are  restricted  in 
form,  but  include  many  of  the  practical  devices  often  encountered  in  detection 
and  estimation  problems.  The  third  class  covers  the  most  general  second-order 
processor;  it  requires  the  solution  for  the  eigenvalue  and  modal  matrices  of 
an  MxM  matrix  (where  M  is  the  size  of  the  general  quadratic  form)  in  addition 
to  the  program  furnished  here.  The  approach  utilized  here  allows  a  user  to 
quickly  and  easily  obtain  accurate  quantitative  information  about  the 
performance  of  a  particular  processor,  and  to  investigate  the  effects  of 
making  changes  in  any  of  the  input  constants  or  parameters. 

Approximations  to  the  performance  of  continuous  quadratic  processors  are 
possible  by  use  of  the  above  procedures.  For  example, 


mAp  nA2)  x(nA2),  (101) 


which  is  of  the  form  X^BX  encountered  in  (62).  Also 


dt^  dt2  u(t^)  e(tp  t2)  v(t2)  2  a^  a2  </_u(mA1)  e( 


mAp  nA?)  v(nA2),  (101) 


m,n 


T  12 

which  is  of  the  form  U  B  V  encountered  in  (79). 
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Receiver  operating  characteristics,  that  is,  detection  probability  vs 
false  alarm  probability,  can  be  easily  determined  from  the  above  results. 
First  store  the  exceedance  distribution  for  zero  signal  strength  in  an  array. 
Then  plot  the  exceedance  distributions  for  nonzero  signal  strengths  vs  this 
stored  array  of  numbers,  each  point  for  a  common  threshold.  The  common 
thresholds  are  most  easily  realized  by  keeping  sampling  increment  a  and  FFT 
size  the  same  throughout  all  the  computations. 
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APPENDIX  A.  SECOND-ORDER  PROCESSOR 


This  program  computes  the  cumulative  and  exceedance  distribution 
functions  of  random  variable  (8)  via  characteristic  function  (9).  The 
required  inputs  are  listed  in  lines  20-120  and  are  annotated  consistently 
with  (8).  The  parameters  Dp  Dp  Ng,  Np  N£  required  in  characteristic 
function  (9)  are  pre-computed  once  in  loop  290-510  for  the  sake  of  execution 
time.  The  mean  of  x  is  entered  in  line  520.  When  we  enter  loop  590-830  for 
the  actual  calculation  of  the  characteristic  function  (9),  the  number  of 
computations  are  minimized.  For  example,  only  one  complex  exponential  and 
square  root  are  required  per  f  value,  in  lines  740-750.  The  square  root  in 
(9)  is  not  a  principal  value  square  root,  but  in  fact  must  yield  a  continuous 
function  in?.  In  order  to  achieve  this,  the  argument  of  the  square  root  is 
traced  continuously  from  0  (line  530).  If  an  abrupt  change  in  phase  is 
detected,  a  polarity  indicator  takes  note  of  this  fact  (line  780)  and  corrects 
the  final  values  of  characteristic  function  f  (f)  (line  790).  More  detail 
on  the  selection  of  L,  a,  b  in  lines  130-150  is  available  in  [1]. 


10  ! 
20 

SECOND-ORDER  PROCESSOR 

K  =  5 

J 

N  u  m  b  e  r  o  f  t  e  r  m  s.  s  u  rn  rn  e  d 

30 

DATA  . 6, -.5, . 4, 

“  •  3  <t  .2 

j 

-=l  <  k  >  w  e  i  g  h  t  i  n  g  s 

40 

DATA  «  9 ,  . 3 , . 7 , - 

.  6 ,  -  .  5 

! 

b  < k )  w e i g h t i ngs 

58 

DATA  - . 6 , - . 3 , 1 , 

1 . 2,  1 . 4 

1 

•:  <  k  )  w  e  i  g  h  t  i  n  g  s 

63 

DATA  .  1 ,  -  . 2, 3 

, .4,  .5 

j 

d <  k )  weight! n g  s 

70 

DATA  7, . 6, .5, 

.  4 ,  - .  3 

| 

e  <  k  )  w  e  1  g  h  t  i  n  g  $ 

8  0 

DATA  . 2, . 3, . 4, - 

.  5  ,  -  .  6 

1 

M  e  an  s  o  f  r  an  d  o  m  y  an  i  ab  1  e  s  s  <.  k  ') 

90 

DATA  .3,-. 7, -.6 

, .5, .4 

I 

M  e  an  s  o  f  r  an  d  o  m  y  ar  l  a  b  1  e  s  t  <  k  ) 

100 

DATA  .  1  ,  . 3, . 5,  . 

7,  .9 

1 

S  t  an  d  ar  d  d  e  y  i  at  i  o  ri  •=•  o  f  s  ( k  ) 

110 

DATA  . 2, . 4, . 6, . 

3,  1 

I 

S  t  an  d  ar  d  d  e  y  i  at  i  o  n  s  o  f  t  k  ) 

120 

130 

140 

DATA  . 4 , — . 5 ,  . 6 , 
L  =  25 

Del  t  a  —  ■  0  5 

.  7  ,  -  .  8 

1 

1 

j 

Cornel  at i o n  c o e f f s .  o f  s C k )  an d  t  < k ) 
Limit  on  integral  o f  c h ar ,  f u n c t i o n 

S  am  p  1  i  n  g  incre  rn  e  n  t  o  n  c  h  ar  .  f  u  n  c  t  i  o  n 

150 

160 

170 

Bs=. 75*<2*PI/De 
rif=2^3 

PRINTER  IS  0 

1 1  a ) 

j 

j 

Shift  b,  as  fraction  of  alias  interval 
Size  of  FFT 

1S0  PRINT  "L  =";L, "Delta  =";Belta,"b  =";Bs,"Mf  =";Mf 

190  REDIM  fta:K),B<l:K),C(l:K),D(i:K),E(l:K) 

200  RED  I M  Ms  < 1 : K  >  ,  Mt  < 1 : K ) , Ss  < 1 : K  > , St  < 1 : K> , Rho< 1 ! K) 

210  REDIM  D 1 < 1 : K  ) , D2  < 1 : K ) , N0  <l:K),Nl<l:K),N2<l:K> 

220  REDIM  X<0:Mf-l),Y<0:Mf-l) 

230  DIM  fl a : 1 0 > , B < l : 1 0 ) , C a : 1 0  ) , D < 1 : 1 0 > , E < 1 : 1 0 > 

240  dim  Ms <  1 : 1 0 ) ,  Mt  <  1 : 1 0 ) ,  Ss <  l :  1 0  :> ,  St  <  1 : 1 0  > ,  Rho <  1 : 1 0  ) 
250  DIM  D 1< 1 :  1  0 > , D2 a : 1 0 ) , N0 <1: 1 0 >, N 1 < 1 : 10),N2<1: 10) 
260  DIM  X<0: 1023) , Y<0: 1023) 

270  READ  fl  <  * ) , E ( *  > , C  <  * > ,  D ( * ) ?  E < * )  !  Enter 

280  R  E  fl  D  M  s  *'.  *■)  j  M  t  (  *  )  ,  S  s  (  *  ) ,  S  t  (  *  )  f  F!  h  o  (  *  )  *  c  o  n  s  t  -an  t  s 
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290  FOR  J= 1  TO  K  !  Calculation 

300  T 1  =Ms  <  J  >  ’-“2  !  of 

310  T  2  =  M  t  <  J ) A 2  !  p a r am  e-  t  e  r  s 

320  T3-Ss  <  J  >  A2 

330  T  4  =  S t < J)  A2 

340  T5-Ms  <  J  >  *Mt < J) 

350  T  6 = R  h  o ( J  >  *  S  s < J ) *  S t ( J ) 

360  T7=4*A  <:  J  >  *B  (  J  >  -C  <  J  >  A2 
370  T  8  =  <1 - R  h  o <  J  >  A2 ) *T3*T4 
3S0  T9=Mt < J)*Ss< J) 

390  T 10  =  Ms<J)*St <  J ) 

400  T 1 1 =D  <  J ) A2 

410  T 12-E(J)A2 

420  T  1  3  =  D  (  J  )  *  E  (  J  > 

430  D1 < J>=2*<fl< J>*T3+B< J)*T4+C< J>*T6> 

440  D2<J>=T7*T3 


450 

460 

470 

430 

490 

500 

510 

529 

530 
540 
550 
560 
570 
580 
590 
600 
610 
620 
630 
640 
650 
660 
670 
638 
690 
700 
710 
720 
730 
740 
750 
760 


N0<  J  )  =fl <  J ) *T  1 +B <  J >*T2+C <  -J  >  *T5  +  B <  J )  *Ms  <  J )+E<  J )  *M t  <  J > 
T  =  T7*  < . 5*<T1*T4+T2*T3)-T5*T6> 

T  =  T+  <  2*A  < J ) *E  <  J  > -C < J  >*D< J > ) *Ss< J  > *  <T9-Rho < J ) *T 1 0 > 

T  =  T+  < 2*B  <  J ) *D  < J  > -C <  J >  *E  <  J  > )*St ( J > *<T1 0-Rho  < J > *T9 > 

HI < J)=T-. 5*<T11*T3+T12*T4>-T13*T6 


N2  < J ) =- <  fl < J>*T 1 2  +  B  <  J  >  *T 11 -C  < J ) *T 1 C 
NEXT  J 

riux=SUM  <  N0  >  +  .  5*3 U M  <  D 1  > 

R  =  0 
P=1 

Muy =Mux+Bs 
X<0)=0 

V  <  0  >  =  .  5  *  D  e  1  t  a* M u y 
N  =  I  N  T  <  L  D  e  1  t  a  ) 

FOR  Ns= 1  TO  N 
X i =  D  e 1 t  a*N s 
X  2  =  X  i  *  X  i 
Pr  =  1 

Pi =  S  r  =  S i =0 
FOR  J=1  TO  K 
D  r  =  1  -  X  2  *  D  2  <  J  ) 

D i =-X i *D 1 ( J ) 

CALL  Mul <Pr, Pi , Dr , Di , A, B > 

Pr=A 
Pi  =B 

CALL  Di u < N 0 < J > - X 2 * N 2 < J > , -Xi *N1 < J) , 
Sr =Sr +A 
Si =S i +  B 
NEXT  J 

CALL  Exp  < -X i *S i  , X i *  <  Sr  +  Bs  > , A , B ) 
CALL  Sqr<Pr,Pi,C,D) 

Ro  =  R 


M  e  an  o  f  r  an  d  o  m  v  ar  i  ab  1  e  x 
Argument  of  square  root 
P  o  1  ar  i  t  y  i  n  d  i  c  at  o r 


A  r  g  u  m  e  n  t  x  i  o  f  c  h  ar .  f  n  . 
Cal  c  u  1  at  i  on 
of 

c  har  ac  t  er  i  s  t,  i  c 
f  u  n  c  t  i  o  n 
f y <xi > 


D  r ,  D  i  ,  A ,  B ; 


770  R  =  ATN  <  D^C  > 

780  IF  ABS(R-Ro)  >1 • 6  THEN  P=-P 
790  CALL  Di y<R, B, C*P, D*P, Fyr, Fyi ) 

800  Ms  =  Ns  MOD  Mf  !  Collapsing 

319  X  <  Ms  > =X < Ms ) +  Fy r/Ns 

320  Y ( Ms ) =Y  <  Ms ) +Fy i /Ns 
830  NEXT  Ns 

340  CALL  Fftl0z<Mf ,X<*>,Y<*>>  !  0  subscript  FFT 
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850 

869 

379 

830 

890 

900 

910 

920 

930 

940 

950 

969 

970 

930 

990 

1000 

1010 

1020 

1  030 

1040 

1050 

1060 

1070 

1030 

1090 

1  100 

1110 

1  120 

1  130 

1  140 

1150 

1  160 

1170 

1  130 

1  190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1230 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 
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PLOTTER  IS  “GRAPHICS" 

GRAPHICS 

SCALE  0„ Mf , -14, 0 
LINE  TYPE  3 
GRID  MfV3,l 
PENUP 

LINE  TYPE  1 

B  =  B  s  *  M  f  *  B  e  1  t  a--'  <  2  * P  I  >  !  0  r  i  q  i  n  f  o  r  r  an d  o m  y  a r  i  ab  1  e  x 

MOVE  B,0 
DRAW  B , - 1 4 
PENUP 

FOR  Ks=0  TO  Mf-1 
T=Y<Ks>/PI-Ks/Mf 

X < Ks >  =  . 5-T  !  Cumu 1  at i ye  probab i 1 i t y  in  X < * > 

Y < K s >  =  P r  = . 5  +  T  !  Exceed an c e  p r oh ab i 1 i t y  in  Y < * ) 

IF  Pr >=1E-12  THEN  Y  =  L G T ( P r > 

IF  PrO- IE- 12  THEN  Y  =  -24-LGT < -Pr > 

IF  ABSCPrX  IE-12  THEN  Y  =  -12 
PLOT  K s , Y 
NEXT  Ks 
PENUP 

PR  I  NT  Y  ( 0 )  ;  Y  (1  >  ;  Y  <  Mf -2  >  ;  Y  Mf-1  > 

FOR  Ks=0  TO  Mf-1 
Pr  =  X  <  Ks ) 

IF  Pr  >= 1 E- 1 2  THEN  Y=LGT(Pr ) 

IF  PrO- IE- 12  THEN  Y  =  -24-LGT < -Pr > 

IF  ABSCPrX  IE- 12  THEN  Y  =  -12 

PLOT  Ks , Y 

NEXT  Ks 

PENUP 

PAUSE 

DUMP  GRAPHICS 
PRINT  L I N < 5 ) 

PRINTER  IS  16 
END 
! 

SUB  M u  1  < X 1 ,  Y 1 ,  X 2 ,  v' 2 ,  A ,  B )  I  Z  1  * Z 2 

A=X 1 *X2-Y 1 *Y2 
B=X 1 *Y2+X2*Y 1 
SUBEND 

i 

SUE  D i y  <  X 1 , Y 1 , X 2 , Y 2 , A , B  >  !  Z1/Z2 

T=X2*X2+Y2*Y2 

A=<X1*X2  +  Y1*Y2XT 

B  =  <Y1*X2-X1*Y2XT 

SUBEND 

i 

SUB  Exp<X,Y,A,B)  !  E X P ( Z ) 

T =EXP  <  X ) 

A=T*CGS<Y) 

B=T  *S I N  <  Y> 

SUBEND 
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1380  SUB  Sqr<X,Y,fl,B>  !  PRINCIPAL  SQR<Z> 

1390  IF  XO0  THEN  1430 
1400  A=B=SQR<. 5*ABS<Y>> 

1410  IF  Y< 0  THEN  B=-B 

1420  GOTO  1540 

1430  F=SQR<SQR<X*X+Y*Y>) 

1440  T=. 5*flTN<Y/X) 

1450  A=F*C03 ( T ) 

1460  B=F*SIN<T) 

1470  IF  X>0  THEN  1540 
1480  T  =  A 
1490  fl=-E 
1500  B=T 

1510  IF  Y >=0  THEN  1540 

1520  fl=-R 

1530  B=-B 

1540  SUBEND 

1550  ! 

1560  SUB  Fftl0z<N,X<*),Y<*>>  !  N  <=  2M0  =  1024,  N=2"- INTEGER  0  subscript 

1570  DIM  C  <  0 : 256) 

1580  INTEGER  II, 12, 13, 14, 15, 16, 17, 18,19, 110, J,K 

1590  DATA  1, .999981175283, .999924701839, .999830581796, .999698818696, .9995294175 

01 , . 999322384588, . 999077727753, . 998795456205, . 998475580573, . 9981181 12900 

1600  DATA  . 997723066644, . 997290456679 , .996820299291, .996312612183, .995767414468 

, . 995184726672, .994564570734, .993906970002, .993211949235, .992479534599 

1610  DATA  . 991709753669, . 990902635428, . 990058210262, . 989176509965, . 988257567731 

,  .  987301418158, . 986308097245, . 985277642389, . 984210092387, . 983105487431 

1620  DATA  . 9819638691 10, . 980735280403, . 979569765685, . 978317370720, . 977028142658 

, .975702130039, .974339382786, .972939952206, .971503890986, .970031253195 

1630  DATA  .968522094274, . 96697647 1 045 , .965394441698, .963776065795, .962121404269 

, . 960430519416, . 953703474896, . 956940335732, . 955141168306, . 953306040354 

1640  DATA  . 951435020969, . 949528180593, . 947585591018, . 94560732538 1 , . 943593458162 

, .941544065183, .939459223602, .937339011913, .935133509939, .932992798335 

1650  DATA  . 930766961079, . 928506080473, . 926210242133, . 923879532511 , . 921514039342 

, . 9191 13851690, . 916679059921 , . 914209755704, . 911706032005, . 909167983091 

1660  DATA  . 906595704515, . 903989293123, . 901348847046, . 898674465694 , . 895966249756 

, . 893224301196, . 890443723245, . 887639620403, . 88479709843 1 , . 831921264348 

1670  DATA  . 879012226429, . 376070094195, . 873094978418, . 870086991 109, . 867046245516 

, . 863972856122, . 860866938638 , . 857728610000, . 354557988365, . 851355193105 

1680  DATA  . 848120344803, . 844853565250, . 841554977437, . 838224705555, . 834862874986 

, . 831469612303, . 828045045258, . 824589302785, . 321 102514991 , . 817584813152 

1690  DATA  .814036329706, . 810457198253, . 806847553544, . 80320753 1 48 1 , . 799537269108 

, . 795836904609, . 792106577300, . 738346427627, . 784556597 1 56 , . 780737228572 

1700  DATA  . 776838465673, . 773010453363, . 769 1 03337646 , . 765 1 67265622 , . 761202385484 

,  .  757208846506, . 753186799044, . 749136394523, . 74505778544 1 , . 740951125355 

1710  DATA  . 736816568377, . 732654271672, . 728464390448, . 72424703295 1 , . 720002507961 

, . 715730825284, .71 1432195745, . 707106781 187, . 702754744457, . 698376249409 

1720  DATA  . 693971460390, .689540544737, . 685033667773 , . 680600997795 , .676092703575 

, . 671553954847, . 666999922304, . 662415777590, . 657806693297, . 653172842954 

1730  DATA  . 648514401022, . 643831542390, . 639124444864, . 634393234 1 64 , . 629638238915 

, . 624859488142, . 62005721 1 763, . 6 1 523 1 59058 1 , . 6 1 0382806276 , . 605511041404 

1740  DATA  . 600616479334, . 595699304492, . 590759701359, . 585797857456, . 530813953096 

, . 575808191418, . 570780745887, . 565731810784, . 56066 1 576 1 97 , . 555570233020 

1750  DATA  . 550457972937, . 545324988422, . 540171472730, . 534997619387, . 529803624636 

, . 524589682678, . 519355990166, .514102744193, . 508830142543, . 503538383726 
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1760  DAT ft  . 498227666973, . 492898192230, .487550160148, .482183772079, . 476799230063 

, . 471396736826, . 465976495768 , . 460538710958, . 455083587126, . 44961 1329655 

1770  HftTfl  . 444122144570, . 4386 1 6238539 , . 433893818853, . 427555093438, . 422808270800 

, . 416429568098, . 410843171858, . 405241314085, . 399624 1 99846 , . 393992840861 

1780  DftTft  . 388345846699, . 382683432365 , . 377007410216, . 371317193952, . 365612997805 

, . 359895036535, . 354163525428, . 348418680249, . 3426607 1 73 1 2 , . 336889853392 

1790  HflTft  . 331 106305768, . 325318292162, . 319502030816, . 313681740399, . 387349648842 

, . 302085949319, . 296158888244, . 298234677254, . 234487537211 , . 278519689385 

1300  DftTft  . 272621355450, . 2667 1 2757475 , . 268794117915, . 254365659685 , . 248927685746 

, . 242930179903, . 237023605994, . 23 1 853 1 8828 1 , . 225033911360, .219101248157 

1310  DftTft  .2131 10319916, . 287111376192, . 281184634842, . 195090322016, . 139068664150 

, . 183039887955, . 177004228412, . 170961333768, . 164913120498, . 153353143334 

1828  DftTft  . 152797185258, . 146738474455, . 140658239333, . 134530783587, . 128498110794 

, . 122418675199, . 1 16318638912, . 118222287294, . 184121633872, . 9381 7 1 403296E- 1 

1830  DftTft  . 9 1 903956497 1 E- 1 , . 357973 1 23444E- 1 , . 7968243797 1 4E- 1 , . 735645635997E- 1 , . 

674439 195637E-1 , . 6 1 3287363822E- 1 , . 55 1 952443497E- 1 , . 490676743274E- 1 

1348  DftTft  . 4  2  9  3  8  2  5  6  9  3  4  9  E - 1 ,  . 3638722294 1 4E- 1 , . 30674303 1 766E- 1 , . 2454 1 22S5229E- 1 ,  . 

1 3  4  0  6  7  2  9  9  8  5  8  E - 1 ,  .  1 227 1 5332S57E- 1 ,  . 6  1 35334649 1 5E-2 , 0 

1850  READ  C<*> 

I860  K  = 1 0  2  4 / H 
1370  FOR  J=0  TO  N/4 
1338  C< J>=C<K*J> 

1398  NEXT  J 
1900  N 1 = N / 4 
1910  N2=N 1+1 
1920  N  3  =  N  2  +  1 
1930  N 4  =  N  1  +  N  3 

1948  Log2n= I  NT  < 1 . 4427*L0G < N >  + .  5> 

1958  FOR  11  =  1  TO  L o g 2 n 
I960  I  2  =  2-^  <  Loq2n~  I  1  > 

1970  13=2*12 

1930  I4=N/I3 

1998  FOR  15=1  TO  12 
2000  16= < 15-1 )  *  I  4  +  1 

2810  IF  I  6  <  =  N  2  THEN  2858 
2820  N6=-C<N4-I6-1 > 

2030  N  7  =  -  C  ( I  6  -  N  1  -  1 ) 

2040  GOTO  2078 
2050  N6=C<I6-1> 

2068  N  7  =  -  C  ( N  3  - 1  6  - 1 ) 

2078  FOR  17=8  TO  H-I3  STEP  13 
2880  13=17+15 

2098  19=18+12 

2100  N3=X<I8-l)-X<I9-i> 

21  10  N9  =  V< 18-1 >-Y< 19-1 ) 

2128  X< 13-1 )=X< I  3-1 >+X<I9-l > 

2130  Y< 13- 1>=Y( 18-1 )+Y (19-1 > 

2140  X ( 19-1 >=N6*N8-N7*N9 
2158  Y  < 19-1 )  =  N  6  *  N  9  +  N  7  *  N  8 
2160  NEXT  17 
2178  NEXT  15 
2138  NEXT  II 
2198  1 1 =Log2n+ 1 

2280  FOR  12=1  TO  10  !  2M  0=1 824 

2210  C< 12—1 )=1 

2220  IF  1 2  >Log2n  THEN  2248 

2230  C< 12-1 >=2A< 11-12) 

2248  NEXT  12 
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2250 

K=  1 

2260 

FOR  11  =  1  TO  C < 9 > 

2270 

FOR  12=11  TO  C < 8 ) 

STEP 

C  <  9  > 

22S0 

FOR  13=12  TO  C < 7 > 

STEP 

C  <  8  > 

2290 

FOR  14=13  TO  C  <  6  > 

STEP 

C<7> 

2300 

FOR  15=14  TO  C < 5 > 

STEP 

C  <  6  > 

2310 

FOR  16=15  TO  C < 4 > 

STEP 

C  <  5  > 

2320 

FOR  17=16  TO  C<3> 

STEP 

C  (  4  ) 

2330 

FOR  13=17  TO  C  <  2 ) 

STEP 

C  ( 3  ) 

2340 

FOR  19=13  TO  C(l) 

STEP 

C  <  2  > 

2350 

FOR  110=19  TO  C<0> 

STEP 

C<1> 

2360 

J=I  10 

2370 

IF  K> J  THEN  2440 

2380 

fl=x<K-n 

2390 

X<K-1 )=X< J- 1 ) 

2400 

X< J-l >=R 

2410 

fl=Y  <K- 1 ) 

2420 

Y  <  K  - 1  >  =  Y  <  J  - 1  > 

2430 

Y< J-l >=R 

2440 

K=K+ 1 

2450 

NEXT  110 

2460 

NEXT  19 

2470 

NEXT  18 

2430 

NEXT  17 

2490 

NEXT  16 

2500 

NEXT  15 

2510 

NEXT  14 

2520 

NEXT  13 

2530 

NEXT  12 

2540 

NEXT  11 

2550 

SUBEND 
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APPENDIX  B.  FADING  FOR  SECOND-ORDER  PROCESSOR 


This  program  computes  the  cumulative  and  exceedance  distribution 
functions  for  characteristic  function  (21),  when  the  power  fading  factor  r  in 

(18)  has  probability  density  function  (19).  The  parameters  Dp  Dp  Nq,  n| 

are  pre-computed  once  in  lines  210-310.  The  logarithms  in  lines  430  and  440 
have  arguments  that  never  cross  the  branch  line  along  the  negative  real  axis 
for  the  principal  value  logarithm;  hence  the  calculated  characteristic 
function  is  automatically  continuous  for  all  ?. 


10  ' 

FADING  FOR 

8  E  C  O  N  D-ORDE  F-I  P  R  0  C  E  S  8  0  R 

20 

N  u  =  2 . 7 

!  Power  1  aw 

for  fading 

30 

K-5 

!  N  u  m  b  e  r  o  f 

terms  summed 

40 

flk  =  .  7 

!  a  <  k  )  w  e  i  g  h  t  i  n  g 

50 

Bk=-. 9 

!  b  <  k  )  w  e  i  g  h  t  i  n  g 

60 

Ck  =- . 6 

!  c < k )  weighting 

70 

DATA  . 2, . 3 

,  .  4 ,  - .  5 ,  - .  6  !  Means  id  f 

r  an  d  o  m  u  an  i  ab  1  e  s  s  <  k  ) 

80 

DATA  .8,-. 

7 , - . 6 , . 5 , . 4  !  Means  o f 

r an  d  o  rn  u  ar  i  ab  1  e  s  t  <  k  ) 

90 

S  s  =  .  3 

!  S  t  an  d  an  d 

deviation  of  s<k) 

1  0  0 

St  =  ,  2 

!  8 1  an  d  an  d 

deviation  of  t<k) 

110 

R  h  o  =  - .  4 

!  Com e  1  at  i  on  c  o e f  f  -  o f  s  <  k  >  an d  t  <  k  > 

120 

L=  1  50 

!  Limit  on 

i  n  t  e  g  r  a 1  o f  c h a r .  f u net i o  n 

130 

Del  t  a= .  25 

!  Samp  ling 

i  n  c  r  e  rn  e  n  t  o  n  c  h  ar  .  f  u  r\  c  t  i  o  n 

140 

Bs= . 625*  <  2 

*  P  I  /  D  e  1 1  a )  !  Shift,  b  , 

as  fraction  of  alias  interv. 

150 

Mf -2*8 

!  Size  of  FFT 

160 

PRINTER  IS 

0 

170 

PRINT  "L  = 

"  ;  L,  “Delta  ="  ;  Delta,  Mb  =,J 

;Bs,"Mf  =M ;Mf 

180 

REDIN  Ms (1 

: K ) , Mt  < 1 : K) , X < 0 : Mf-1  > , Y<0 

: Mf-1) 

190 

DIM  Ms  < 1 :  1 0  > , Mt  <1:1 8  > , X  < 0 : 1 023 )  ,  Y  <  0 

: l 023) 

200 

READ  M s  ( *  ) 

,Mt<*>  ! 

E  n  t  e  r  const  an  t.  s 

210 

M  2  0  =  D  O  T (Ms 

,  Ms)  ! 

Cal  c u  1  at  i  on 

220 

M  0  2  =  D  0  T  <  M  t 

,  Mt)  ! 

of 

230 

Ml  1 =  D  0  T (Ms 

,  Mt)  ! 

p  ar  am  e  t  e  r  s 

240 

T 1 =Ss*Ss 

250 

T2=St*St 

260 

T3=Rho*Ss*St 

270 

T4=4*Ak*Bk 

-Ck  *Ck 

280 

N0p=Ak*M20+Bk*H02+Ck*Ml 1 

290 

N 1 p  =  T4*  < . 5 

*  <  T  2  *  M2  8  +  T  1  *  f*1 0  2  )  -  T  3  *  M  1  1  ) 

300 

Dl=2*<Ak*Tl+Bk*T2+Ck*T3> 

310 

D2=T4*< l-Rho*Rho)*Tl*T2 

320 

D 1  p  =  D  1 +N0p/Nu 

330 

D2p=D2+Nlp/Nu 

340 

M  u  x  =  N  0  p  + .  5  *  K  *  D  1  ! 

Mean  of  random  variable  x 

350 

Muy =Mux+Bs 

360 

T=Nu-. 5*K 

B-l 
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370 

330 

390 

400 

410 

420 

430 

440 

450 

460 

470 

4S0 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

680 

618 

620 

638 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

748 

750 

760 

770 

780 

790 

800 

810 

820 

838 

840 

850 

860 

870 

380 


X  <  0  >  =  0 

Y  <  0  >  = .  5  *  D  e  1 1  a*  M  u  y 
M=  I  NT  <  L/'Be  1 1  a) 

FOR  Ns*l  TO  N 

X i =De 1 t a*Ns 
X  2  =  X  i  *  X  i 

CALL  Log < 1 -X2*B2 , -X i *B 1 , A ,  B) 

CALL  Log < 1 -X2*B2p , -X i *D 1 p , C , D > 

T 1 =T*A-Nu*C 

T2=T*B-Nu*D+Bs*Xi 

CALL  Exp<Tl,T2,Fyr,Fyi) 

M s = N  s  MOB  Mf 
X  <  Ms  )  =X  (  Ms  )  +  Fyr/*Ns 

Y  C  Ms ) =Y  <  Ms  >  +Fy i /Ns 
NEXT  Ns 

CALL  Fft,  10z<Mf  ,X<*>  ,Y(0) 

PLOTTER  IS  "GRAPHICS" 

GRAPHICS 

SCALE  8, Mf , -14, 0 
LINE  TYPE  3 
GRIB  Mf/8,1 
PENUP 

LINE  TYPE  1 
B=Bs*Mf *Bel ta/<2*PI > 

MOVE  B,0 
BRHW  B , -  1  4 
PENUP 

FOR  Ks=8  TO  Mf-1 
T  =  Y  ( K  s  >  /  P  I  -  K  s  /  M  f 
X<Ks>=. 5-T 
Y<Ks)=Pr=. 5+T 

IF  Pr  >= 1 E- 1 2  THEN  Y  =  L  G  T ( P  r ) 

IF  Pr<  =- 1 E- 1 2  THEN  Y  =  -24-LGT < -Pr > 
IF  ABSCPrK  IE-12  THEN  Y  =  -12 
PLOT  K s , Y 
NEXT  Ks 
PENUP 

PR  I  NT  Y  < 0 ) ; Y  < 1 > ; Y  < Mf -2 ) ; Y  < Mf - 1 > 

FOR  Ks=8  TO  Mf-1 

Pr=X<Ks) 

IF  Pr  >= 1 E- 1 2  THEN  Y  =  LGT<Pr> 

IF  Pr<  =- 1 E- 1 2  THEN  Y  =  -24-LGT < -Pr ) 

IF  ABSCF'rX  IE- 12  THEN  Y  =  -12 

PLOT  Ks,Y 

NEXT  Ks 

PENUP 

PAUSE 

DUMP  GRAPHICS 
PRINT  L I N  <  5  > 

PRINTER  IS  16 
END 

I 


!  A  r  g  u  m  e  n  t  x  i  o  f  c  h  a.  r  .  f  n , 
!  Calc  u 1  at  i  o  n 

!  o  f 

!  characteristic 
!  function 
!  f  y  ( x  i  ) 

!  Collapsing 


!  0  subscript  FFT 


0  r  i  g  i  n  f  o  r  r  an d  o m  v  ar  i  ab  1  e  x 


C  u  m  u  1  at  i  u  e  p  r  o  b  ab  i  1  i  t  y  i  n  X  ( 
Exceed an c e  p r o b ab i 1  it y  in  Y < 


B-2 


890 

990 

910 

920 

930 

948 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 


TR  7035 


SUE  ExpCX, Y, R , B)  !  EXF’CZ) 

T=EXP<X> 

A  =  T  *  C  0  S  <  Y  > 

B  =  T  *  S  I N  <  Y  ) 

SUBEND 

j 

SUE  L  o  g  <  X ,  Y ,  ft ,  E  )  !  PRINCIPAL  LOGCZ) 

R  = . 5*LGG<X*X+Y*Y) 

IF  XO0  THEN  1000 
B= . 5*P I *SGH  <  Y ) 

GOTO  1020 
B  =  flTN<  Y/X  > 

IF  X<8  THEN  B  =  B+P I *  < 1-2*  <  Y<  0 ) > 

SUBEND 

I 

SUB  Fft 10z<N , X<*) , Y  <*))  !  N  <=  2--10  =  1824,  N=2A INTEGER 


8  s  u  b  s  c  r  i  p  t 


B-3/B-4 
Reverse  Blank 
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APPENDIX  C.  NARROWBAND  CROSS-CORRELATOR 


This  program  computes  the  cumulative  and  exceedance  distribution 
functions  of  random  variable  (30)  via  characteristic  function  (35).  The 
parameters  D^?  d2,  Nq,  are  pre-computed  in  lines  280-390  and 
weighted  according  to  (35)— (36)  in  lines  400-440.  All  the  functions  employed 
are  analytic. 


10  ! 
20 
30 
40 
50 
60 
70 
80 
90 
100 
1  10 
120 
130 
140 
150 
160 
170 
130 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 

349 

350 
360 
370 
330 
390 
408 
410 
420 
430 
440 
450 
460 


HfiRROWBfiND  CROSS-CORRELATOR 
K  =  5  ! 

D A  T  A  .6, -.5, .4, -.3, .2  ! 

DATA  .9, .3, .7, -.6, -.5  ! 

DATA  -.6, -.3, 1,1.2, 1.4  ! 
DATA  . l,-.2,-.3, .4, .5  ! 

DATA  7, . 6, . 5, . 4, 3  ! 

DATA  .  1 ,  .  3, . 5, .  7,  .  9  ! 

DATA  . 2, . 4, . 6, . 3, 1  ! 

DATA  . 4, 5, . 6, . 7j 8  ! 

DATA  . 9,-. 7,-. 5, . 3, 1  ' 

L=50  ! 

De 1 1  a= ■ 5  ! 

E  s  = . 5  *( 2  *  P I s Del t  a )  ! 

Mf =2*8  ! 


N u m her  o f  t  e  r m s  s u rn m e d 
ui  ( k  >  w  e  i  g  h  t  i  n  g  s 

a 1 < k  >  signal  1  i n - p  h a s e  c o m p o n e n  t  s 

b  1  <  k  >  signal  1  q u ad r at u  r e  c o m p  o n  e n t 

a 2 < k )  signal  2  in- p h  as e  c  o m p o  n e  n  t  s 

b  2  (.  k  >  signal  2  q  u  a  d  r  a  t  u  r  e  c  o  m  p  o  n  e  n  t 

si  g  rn  a  1  <  k  >  noise  1  s  t  an  d  a  r  d  d  e  v  i  at  i  o  n  s 
s  i  g  rn  a2  <  k  >  n  o  i  se  2  s  t  an  id  ar  d  d  e  u  i  at  i  o  n  s 
r h o<k)  noise  i n - p h as e  cor r .  c o e f f s . 

1  am b d a  n o i se  q u ad r at u r e  c o r r .  coe f f s 

Limit  on  i n t e g r a 1  o f  c h ar .  f u n c t i o n 
S  am  p  1  i  n  g  incne  rn  e  n  t  o  n  c  h  ar  .  f  u  n  c  t  i  o  n 
S h i f t  b ,  as  f  ract i on  o f  alias  i n t  e r  o a 1 
Size  of  FFT 


PRINTER  IS  0 


PRINT 
RED  in 
REDIN 
REDIN 
REDIN 
D  I M  W  <  1  : 

DIM  81(1 
DIM  DKl 
DIN  X< 0:1023) 


"L  =" ; L, “Delta 
W<i:K),fli<i:K), 
S l < l : K> 

Dl(l:K) 

x<0: Mf- 


=“ ; Delta, 

B  1  ( 1  :  K  )  ,  A  2  <  1  ; 


'Nf  =  "  ;Mf 


K  > , E2< 1 : K) 


,  S  2  <  1 :  K  > ,  R  h  o  <  1 :  K  > ,  L  am  b  d  a  <  1 :  K  > 

, D2< l : K> , N  0 ( i :  K) , n l < l : K) 


1  > , 


1  0  >  ,  A  1 
:  10) , 

:  10>,  ds 


V<0: Mf-l ) , 
: 1 0 > , B 1 <1 : 


W2 ( l : K) 
10) , A2< 1 


10) , B2< 1 : 10) 


2 <1 : 10) , Rho  < 1:10), Lambda* 1:10) 

<  l :  10) , N0<  l :  10) , Nl< 1 : 10) 

i  <0: 1023) , u 2  < l : io) 


A  1  <  *  )  ,  B  1  (  *  )  , 


A2  *  *  ) 


Rho <*)  ,  Larnbd 


READ  WC*>, 

READ  S  1<  *  )  , S 2  < *  ) 

FOR  J=1  TO  K 
S  1  s  =  S  1  <  J  )  *  2 
S2s=S2< J)*2 
T  1  =3  1  <  J)*S2<:  J) 

DlU)=T2  =  TURho(J) 

D  2  <  J )  = . 25*8 1 s*S2s* < 1 -Rho *  J) 
T3=A1< J)*A2< J)+B1 < J)*B2( J) 


, B2<*) 
<*) 


E  n  t  e  r 
c  o  n  s  t  an  t.  s 
Cal c  u 1  at i on 
of 

par  arn  eters 


S2-L am  b  d a  < J ) A 2 ) 


N  0  <  J )  = . 5*T3 

T 4  =  A2  <  J ) *B  1  <  J)-AK  J)*B2<  J) 

T  5  =  S  2  s  *  <  A  1  <J)*2+B1  <  J  )  -s-2 )  +S  1  s*  <  A  2  <  J  )  -"•2+B2  <  J  )  -"*2 ) 
N 1 <  J )  =  . 1 2 5  *  < T  5 - 2  *  T  2  *  T  3 - 2  *  T 1  *  L am bda<J)*T4) 

NEXT  -J 
NAT  W 2  =  W .  W 
NAT  B 1 =  W . D 1 
NAT  D 2  =  W 2 . D 2 
NAT  N 0  =  W . N 0 
NAT  N1=W2.N1 


N u x  —  S U M CN0^+SUM *  D  1  )  !  M e a n  o f  r an d o rn  o ar  i  a b  1  e  o 

Nuy =Mux+Bs 


Ol 
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47Q  X < 0 >  =0 

4S0  Y  < 0 )  =  . 5*Del ta*Muy 

490  N  =  I N  T  <  L  /  D  e  1 1 a ) 

500  FOR  N  s  =  1  TO  N 

510  Xi  =  Bel  t, a*Ns  ! 

520  X2=Xi*Xi  ! 

530  Pr= 1  i 

540  Pi =Sr=Si =0  ! 

550  FOR  J=1  TO  K  ! 

560  Dr=l+X2*D2< J)  ! 

570  Di =-X i *B1 < J> 

530  CRLL  Mu  1  <  Pr , P i , Dr , D i , ft , B  > 

590  Pr=fi 

600  Pi =B 

610  C ft L L  Diu<N0<J>JXi*NllvJ>,Dr,DiJfl,B> 

620  S  r = S  r  +  ft 

630  Si =  S  i +  B 

640  NEXT  J 

650  CALL  Exp<-Xi *Si , Xi *(Sr+Es) ,H,B) 

660  CALL  Diu<A,B,Pr,Pi,Fyr,Fyi> 

670  M s  =  Hs  MOD  Mf  .  ! 

680  X  <  Ms  > =X  < Ms  > +Fyr/Ns 

690  Y < Ms )  =  Y < Ms  >  +Fy i /Ms 

700  NEXT  Ns 

710  CALL  Fft,  10z<Mf,X<*>,  Y<*>)  ! 

720  PLOTTER  IS  “GRAPHICS" 

730  GRAPHICS 

740  SCALE  0,Mf,-14,0 

750  LINE  TYPE  3 

760  GRID  Mf/8,1 

770  PENUP 

730  LINE  TYPE  1 

790  B=Bs*Mf*Delta/<2*PI)  ! 

800  MOVE  B,0 

810  BRAN  B  ,  -  1  4 

820  PENUP 

830  FOR  Ks=0  TO  Mf-1 

840  T  =  Y  <  K s ) x P I - K s x M f 

850  X  <  K  s  )  =  . 5 - T  ! 

360  Y  <  Ks )  =Pr  =  •  5  +  T  ! 

870  IF  Pr >= 1 E- 1 2  THEN  Y  =  LGT <Pr  ) 

380  IF  Pr<=-1E-12  THEN  Y=-24-LGT < -Pr > 

390  IF  RBS  <  Pr  X 1 E- 1 2  THEN  Y  =  -12 

900  PLOT  Ks,Y 

910  NEXT  Ks 

920  PENUP 

930  PRINT  Y  < 0 ) ; Y  < 1 > ; Y ( Mf -2  > ; Y  <  Mf-1 > 

940  FOR  Ks=0  TO  Mf-1 

950  Pr  =  X  <  Ks ) 

968  IF  Pr>= IE- 12  THEN  Y = L G T < P r > 

970  IF  Pr<  =- 1 E- 1 2  THEN  Y=-24-LGT < -Pr ) 

980  IF  ABS<Pr  X  IE- 12  THEN  Y  =  -12 

990  PLOT  Ks ,  Y 

1000  NEXT  Ks 

1010  PENUP 

1020  PAUSE 


Argument  xi  of  char.  fn. 
Cal  c  u  1  at  i  on 
of 

character i st i c 
funct i on 

f  y  <  x  i  ') 


Col  1  aps  i n g 


0  subscript  FFT 


Or i g i n  f or  random  var i ab  1 e  v 


C u m u 1  at i u e  p r o b ab i 1  it y  i n  X C  *  ) 
E  x c  e  e  d an  c e  p r  ob  ab i 1  it  y  i n  Y ( *  ) 


C-2 


1030 

1040 

1050 

1060 

1070 

1030 

1090 

1100 

1110 

1120 

1130 

1140 

1  150 

1168 

1170 

1130 

1190 

1200 

1210 

1220 

1230 

1240 

1250 
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DUMP  GRAPHICS 
PRINT  LINO) 

PRINTER  IS  16 
END 

j 

SUE  Mul <X1, Y1,X2, Y2,fi,B>  !  Z1*Z2 

A  =  X 1 *X2-Y 1 *Y2 
E  =  X 1  *  Y  2  +  X  2  *  V  1 
SUBEND 
! 

SUE  DiyCXl, Y1,X2,Y2,R,B>  !  Z1/Z2 

T=X2*X2+Y2*Y2 
A=<X1*X2+Y1*Y2>/T 
B=<Y1*X2-X1*Y2>/T 
SUBEND 
! 

SUB  Exp  <  X , Y , fi , B  >  !  EXP<Z> 

T  =  E  X  P  <  X  ) 
fl=T*C0S< Y) 

B=T*S I N  <  Y  > 

SIJEEND 

i 

SUB  Fft 10z<N,X<*),Y<*>>  !  N  <=  2M0  =  1024,  N=2* INTEGER  0  subscript 


C-3/C-4 
Reverse  Blank 
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APPENDIX  D.  REDUCED  QUADRATIC  AND  LINEAR  FORM 


This  program  computes  the  cumulative  and  exceedance  distribution 
functions  of  random  variable  (65)  via  characteristic  function  (66).  The 
required  inputs  to  the  program  are  M  and  the  [xm]f  fdm],  of  (68). 

The  square  root  in  (66)  must  again  be  continuous  and  is  handled  exactly  as  in 
appendix  A.  The  parameters  required  in  the  exponential  of  (66)  are 
pre-computed  in  lines  170-210,  and  the  mean  of  q  is  entered  in  line  220. 


10  ! 
20 
30 
40 
50 
60 
70 
30 
90 
100 
110 
120 
130 
140 
150 
160 
170 
130 
1 90 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
338 
340 
350 
360 
378 
380 
390 
400 
410 


REDUCED  QURDRRT I C  RND  LI  HERR  FORM 


M  =  5 

DATA  .  2  ,  - .  3 ,  .  4 ,  .  5 ,  - .  6 
BfiTR  - . 1, .3, .5, . ? , - . 9 
DRTR  .  6, . 5, 4, 3, .2 
L  =  8  0  0 
Delta=. 03 

Es=.  5*<2*PI/Del  ta> 

M  f  =  2  M  0 


Number  of  terms  summe 
Lambda  val ues 
d  Mai ues 
N  u  values 

Limit  on  i n t  e g r a 1  o f 
S  am  p  1  i  n  g  i  n  c  r  e  m  e  n  t  o  n 
Sh i ft  b ,  as  f rac t i on 
Size  of  FFT 


d 


c  h  ar .  f  u  n  c  t  i  o  n 
char,  function 
of  alias  interv 


al 


PRINTER  IS  0 

PRINT  11 L  ="  ;L,  "Delta  =H;  Delta,  "b  =“  ;Bs,"Mf  =  "  ;  Mf 
RED  I M  Lambda< 1 :  M  > ,  D  <  1 :  M ) , Nu < 1 : M ) , R  < 1 : M > , B < 1 : M > , C  < l : M ) 

RED  I M  X < 0 : Mf - 1 > , V  < 0 : Mf - 1 ) 

D  I  M  Lambda-;  l :  1 0  >  ,  D  <1 :  1 0  >  ,  Nu  <1 : 1 0  > ,  R  a  :  1 0  )  ,  B  <  1 : 1  0  >  ,  C  <  1 :  1 0  ) 
DIM  X<0:  1023) , Y <  0 : 1023) 

R E R D  L  am  b  d  a  <  * ) , D  <  *  > ,  N  u  <  *  >  !  Enter  constants 

FOR  Ms=l  TO  M  !  Calculation 

fl<Ms)=2*Lambda< M s )  !  of  parameters 

B  <  Ms  )  =  <  Larn b d a <  Ms )  *Nu  <  Ms )  +D  <  Ms  >  )  *Nu  (Ms) 

C(Ms)=.5*D(Ms)A2 
NEXT  Ms 


M  u  q  =  S  U  M  <  L  am  b  d  a )  +  S  U  M  <  B ) 

R  =  0 
P=1 

Muy =Muq+Bs 
X(0)=0 

Y  ( 0 ) - .  5  *  D  e 1 1  a*M u y 
N= I  NT ( L/De 1 t  a) 

FOR  Ns= 1  TO  N 
X i =De 1 t a* Ns 
Pr  =  1 

P i =Sr=S i =0 
FOR  Ms= 1  TO  M 
T =-R  <  Ms ) *X i 

CALL  Mul <Pr, Pi , 1 , T, fl, B) 

Pr  =  fl 
Pi  =B 

CRLL  D i m ( B ( M s ) , C(Ms)*Xi , 1 , T , R , B ) 
Sr*Sr+fl 
b  i  —  b  i  +  B 
NEXT  Ms 


Mean  of  random  variable  q 
Argument  of  square  root 
P  o 1 ar it  y  i  n  d i c at  o r 


Argument  x i  of  char.  fn. 
Cal c  u 1  at i on 
of 

c  h ar  ac  t  e r  i  st  i  c 
f  u  n  c  t  i  o  n 
f  y  (  x  i  ) 


D-l 
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420  CALL  E x p  (  -  S  i  *  X  i , ( S r  +  E s ) * X i , A , E ) 

430  CALL  Sqr <Pr , Pi , C, D) 

44Q  Ro=R 

450  R  =  A  T  N  <  D  x  C ) 

460  IF  AES  <  R-Ro ) >1.6  THEN  P  =  -P 
470  CALL  Diy<A,E,C*P,D*P,Fyr,Fyi ) 

430  Ms=Ns  MOD  Mf  !  Collapsing 

490  X  <  Ms  )  =X  <  Ms  )  +Fyr/'Ns  "  j 

500  Y  <  Ms ) = Y ( Ms ) +Fy i xNs 
510  NEXT  Ns 

520  CALL  F  f  t  1  0  z  <  M  f ,  X(*),Y(*))  !  0  subscript  FFT 

530  PLOTTER  IS  "GRAPHICS" 

540  GRAPHICS 

550  SCALE  0 , M f , - 1 4 , 0 

560  LINE  TYPE  3 

570  GRID  Mf/8,1 

530  PENUP 

590  LINE  TYPE  1 

600  E  =  E  s  * M  f  *  D  e 1 1  ax  ( 2  *  P  I )  !  0 r i q i n  f  o r  r an  d  o m  y ar i ab 1 e  q 

610  MOVE  B , 0 

620  DRAW  B , - 1 4 

630  PENUP 

640  FOR  Ks=0  TO  Mf-1 

650  T=YCKs)/PI-Ks/Mf 

660  X<Ks)=.5-T  ’  !  Cum u  1  at  i  y ■=•  p r o b ab  i  1  i  t,  y 

670  Y  <  K  s  )  =  P  r  =  .  5  +  T  !  E  x c  e  e  d  an  c  e  p  r  o  b  ab  i  1  i  t  y 

680  IF  Pr  >= 1 E- 1 2  THEN  Y  =  LGT < Pr ) 

690  IF  Pr <  =- 1 E- 1 2  THEN  Y=-24-LGT < -Pr > 

700  IF  ABSCPrX  IE-12  THEN  Y  =  -12 
710  PLOT  Ks,Y 

720  NEXT  Ks 

730  PENUP 

740  PRINT  Y  < 0 ) ; Y (1 > ; Y  < Mf -2  > ; Y  < Mf-1 > 

750  FOR  Ks=0  TO  Mf-1 
760  Pr=X<Ks) 

770  IF  Pr>=lE-12  THEN  Y  =  LGT <Pr ) 

780  IF  PrO-lE-12  THEN  Y=-24-LGT < -Pr > 

790  IF  ABS(Pr>< IE-12  THEN  Y=-12 
800  PLOT  Ks , Y 

810  NEXT  Ks 

820  PENUP 

830  PAUSE 

840  DUMP  GRAPHICS 

850  PRINT  LINO) 

860  PRINTER  IS  16 

870  END 

830  ! 


i n  X 
i  n  Y  <  *  > 
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890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1030 

1090 

1  100 

1  110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1  190 

1200 

1210 

1220 

1230 

1240 


TR  7035 


SUB  Mul <X1, Y1,X2,Y2, A,B> 
A=X1*X2-Y1*Y2 

B=X1*Y2+X2*Y1 

SUBEND 

i 

!  Z 1 *Z2 

SUB  Di  mOU  ,  Y1  ,X2,  Y2,  ft,  E) 

T  =  X  2  *  X  2  +  Y  2  *  Y  2 

fl=<Xl*X2+Yl*Y2)/T 

B=<Y1*X2-X1*Y2>^T 

SUBEND 

!  Z1/Z2 

SUB  Exp(X,Y,A,B> 

T  =  EXP <  X  ) 

A  =  T  *  C  G  S  (  Y  ) 

B=T*S I N  <  Y ) 

SUBEND 

i 

rj 

Ol 

X 

LU 

SUB  S q r  <  X , Y , A , B  > 

IF  XO0  THEN  1110 

A  =  E  =  SQR  < ■ 5*ABS<Y>> 

IF  Y<  0  THEN  E  =  -B 

GOTO  1220 

F=SQR<SQR<X*X+Y*Y> ) 

T=. 5*ATN<Y/X> 

A=F*COS ( T ) 

B  =  F*S I N  <  T ) 

IF  X>0  THEN  1220 

T  =  A 

A  =  -B 

B  =  T 

IF  Y  >  =  0  THEN  1220 

A*-fl 

B  =  -B 

SUBEND 

i 

!  PRINCIPAL  SQR <  Z > 

SUB  Fft 10z<N,X<*>,Y<*>>  !  N  <= 

2 A  10  =  1024,  H  =  2 A  I  N  T  E  G  E  R  0  s  u  b  s  ■:  r  i  p  t. 
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230 

FOR  1=1  TO  Tt 

290 

A  =  0 

300 

FOR  J=1  TO  K 

310 

V 1  =  R  N  D  -  .  5  ! 

GENERATE  TWO 

320 

V2=RND- . 5  ! 

INDEPENDENT 

330 

S=V 1 *V 1 +V2*V2  ! 

GAUSSIAN 

340 

IF  S  >  • 25  THEN  310  ! 

RANDOM 

350 

Q=<L-LOG<S> ! 

VARIABLES  VIA 

360 

Q  =  3QR  <  Q  +  Q )  ! 

ACCEPTANCE 

370  . 

G1=V1*Q  ! 

AND 

380 

G2=V2*Q  ! 

REJECTION 

390 

S=Ms <  J ) +  A 1 < J  > *G 1 +Ee < J > *G2 

400 

T  =  Mt  < J)+St < J)*G1 

410 

X-X+FK  J>*S*S+BCJ>*T*T+C<J)*S*T  +  D<J)*S+E<:J>*T 

420 

NEXT  J 

430 

xa  >=x 

440 

NEXT  I 

450 

MAT  SORT  X 

460 

PLOTTER  IS  "GRAPHICS'* 

470 

GRAPHICS 

480 

SCALE  -30, 10, -4,0 

490 

GRID  5, 1 

500 

PENUP 

510 

FOR  1=1  TO  Tt 

520 

Y  =  LGT<  < I-. 5>/Tt > 

530 

PLOT  X  <  I  )  ,  Y 

540 

NEXT  I 

550 

PENUP 

560 

FOR  1=1  TO  Tt 

570 

Y=LGT< 1-<I-. 5>/Tt ) 

580 

PLOT  X < I > , Y 

590 

NEXT  I 

600 

PENUP 

610 

END 

E-2 
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